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Abstract 

We consider a Hamiltonian system with 2 degrees of freedom, with a hyperbolic equilibrium point having a loop 
or homoclinic orbit (or, alternatively, two hyperbolic equilibrium points connected by a heteroclinic orbit), as a 
step towards understanding the behavior of nearly-integrable Hamiltonians near double resonances. We provide a 
constructive approach to study whether the unstable and stable invariant manifolds of the hyperbolic point intersect 
transversely along the loop, inside their common energy level. For the system considered, we establish a necessary 
and sufficient condition for the transversality, in terms of a Riccati equation whose solutions give the slope of the 
invariant manifolds in a direction transverse to the loop. The key point of our approach is to write the invariant 
manifolds in terms of generating functions, which are solutions of the Hamilton-Jacobi equation. In some examples, 
we show that it is enough to analyse the phase portrait of the Riccati equation without solving it explicitly. Finally, 
we consider an analogous problem in a perturbative situation. If the invariant manifolds of the unperturbed loop 
coincide, we have a problem of splitting of separatrices. In this case, the Riccati equation is replaced by a Mel'nikov 
potential defined as an integral, providing a condition for the existence of a perturbed loop and its transversality. 
This is also illustrated with a concrete example. 

Keywords: transverse homoclinic orbits, hyperbolic equilibria, Hamilton-Jacobi equation, Riccati equations, split- 
ting of separatrices, Mel'nikov integrals. 



1 Introduction 



1.1 Setup and main results 

The study of the behavior of a Hamiltonian system near a double resonance is one of the main difficulties related 
with Arnol'd diffusion, a phenomenon of instability in perturbations of integrable Hamiltonian systems with more 
than 2 degrees of freedom. Such a behavior is usually studied with the help of resonant normal forms. Neglecting 
the remainder, the normal form can be reduced to a Hamiltonian with 2 degrees of freedom, that in general is 
not integrable. As a first step towards studying the complete system near the resonance, a good understanding of 
this reduced Hamiltonian is very important, and particularly the intersections between the invariant manifolds of 
equilibrium points, along a homoclinic orbit. 

"The authors were supported in part by the Spanish MICINN-FEDER grant MTM2009-06973, and the Catalan CUR-DIUE grant 
2009SGR859. Besides, the author JRP was also supported by MICINN-PN (I+D+I) 2008-2011 (grant JC2009-00306). 
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1.2 



As a model for the reduced system, we consider a classical Hamiltonian with 2 degrees of freedom, of the type 
kinetic energy plus potential energy, with the 2-dimensional torus T 2 as the configuration manifold (see Section 
for more details). In fact, our approach will also be valid in a more general manifold. 

Let H be a Hamiltonian with 2 degrees of freedom, defined on a phase space T*Q, where Q is a 2-dimensional 
configuration manifold. Let q = (91,(72) £ C C R 2 be some local coordinates for Q. Then, we have canonical 
coordinates x = (q,p) € U xl 2 for T*Q, with the standard symplectic form fl = dq A dp, whose associated matrix is 

J = ( ^ , ^ n I . In these coordinates, our Hamiltonian takes the form 



-Id 



H(q,p) = UB(q)p,p)+V(q), (1) 



with a positive definite (symmetric) matrix function B{q), and a scalar function V(q), providing the kinetic energy 
and the potential respectively. The Hamiltonian equations are x = Xjj(x) = JVH(x) = ( ^p, — 4^ J , namely 



I dp ' dq 

q = B(q)p, p = -^ q [{B(q)p,p)]-VV(q), (2) 

We assume that B(q) and V(q) are smooth functions on U , i.e. they are C r with r > 2, or analytic. Then, the 
Hamiltonian equations ^ are C r ~ 1 or analytic. 

For a given homoclinic orbit or loop, biasymptotic to a hyperbolic equilibrium point, our goal is to provide a 
constructive approach to study the transversality of the invariant manifolds along that orbit, inside the energy level 
where they are contained. In fact, we develop our aproach for the case of a heteroclinic orbit, which makes no difference 
with respect to the homoclinic case. We denote O, O two (possibly equal) hyperbolic equilibrium points, and W u ' s , 
W u,s their respective unstable and stable invariant manifolds. Let 7 be a heteroclinic (or homoclinic) orbit, that we 
assume known, connecting the two points, i.e. 7 C W u f~l W s , and we we have to study whether such intersection is 
transverse. 

We consider an open neighborhood U of the first hyperbolic point O, with coordinates (q,p) E U x W 2 as in ([I]) . Of 
course, this neighborhood U may not contain the whole orbit 7, nor the second point O. For the point O, we consider 
a neighborhood U with coordinates (q,p) <E U x M 2 . We assume that the neighborhoods U and U have intersection, 
in which the symplectic change between the coordinates (q,p) and (q,p) is induced by a change in the configuration 
manifold Q: 

q = x(q), P = Vx(q)~ T p (3) 

(where T>x(q) stands for the Jacobian matrix of the change, and we use the notation A~ T for the inverse of the 
transpose of a matrix A). In the intersection U C\U, we will study the transversality between the unstable manifold 
W u of O and the stable manifold W s of O along the orbit 7. 

The transversality between the invariant manifolds will be studied in the coordinates (q,p) oiU. When restricted 
to the neighborhood U, we may refer to the 'outgoing parts' of 7 and W u as the local outgoing orbit and the local 
unstable manifold respectively, before leaving (forward in time) the neighborhood U. On the other hand, since the 
global manifold W s contains the whole orbit 7, it also enters (backward in time) in the neighborhood U, and will be 
compared with W u . 

Let us describe our hypotheses of the Hamiltonian ([I]), expressed in the coordinates (q,p) of the neighborhood U. 
First, we assume: 

(HI) the potential V(q) has a nondegenerate maximum at q = (0, 0), with V(0, 0) = 0. 

Thus, we are assuming that the point O is the origin (q,p) = (0, 0, 0, 0), and this hypothesis says that O is a hyperbolic 
equilibrium point of the Hamiltonian H . Notice that the orbit 7 is then contained in the zero energy level of H . For 
the sake of simplicity, we also assume: 

(H2) the outgoing part of 7 satisfies (72 = 0, with qi increasing along the orbit; 
(H3) the expansion of B(q±,q2) in qi has no term of order one, i.e. f^(<7i,0) = 0. 
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In fact, there always exist local coordinates q — (91,92) such that (H2) is satisfied, though it may be difficult 
to construct them explicitly in a concrete case (nevertheless, see the example of Section |4.2| ). We point out that 
hypothesis (H2) imposes strong restrictions on the form of the functions B(q) and V(q) in ([!]), that define the 



Hamiltonian H in the chosen coordinates (q,p) (see Section 2.3). Concerning (H3), this is not an essential hypothesis, 



but we assume it just to ease the computations, for it is satisfied in all the examples considered in this paper (see 
technical remarks in Sections 12.31 and 12.41) . 



Analogous hypotheses to (HI H3) could be formulated in the local coordinates q = (qi, q 2 ), concerning the second 
point O and the ingoing part of 7. Notice that imposing hypothesis (H2) for both the coordinates q and q implies 
that the change Q satisfies the relation x(<li: 0) = (Qi > 0)- We will express this relation as Xo(<7i) = <7i- Nevertheless, if 
some reversibility properties can be applied, it will be enough to impose hypotheses (HI H3) only for the coordinates 
q, since they imply the analogous ones for q. 

To establish the transversality along 7, we shall express the invariant manifolds in terms of generating funcions. We 
show in Section [2~T] that , in a neighborhood of the point O, the unstable manifold W u can be seen, in the coordinates 
(q,p), as a graph of the form 

p = VS u (q), qeU. (4) 

This gradient form is closely related to the Lagrangian properties of the invariant manifolds. The generating function 
S u (q) can be extended along a neighborhood of the orbit 7 as far as the unstable manifold can be expressed as a graph. 
We will assume that the neighborhood U is such that the form Q is valid for the whole part of W u inside U. The 
same can be done for the stable manifold W s of the second point O. In the coordinates (q,p) of the neighborhood U, 
this manifold also becomes a graph 

p = VS s (q), qeU. (5) 

We shall assume that the neighborhoods intersect, and both generating functions can be extended up to this intersec- 
tion: 

(H4) there exists a point P e 7 <~)U DU with ^-coordinates (q*,Q) £ U and (q*,0) = x(?*>0) 6 U, such that the 
generating functions in Q and ^ can be extended along the outgoing and ingoing parts of 7, up to q* and q\ 
respectively. 

Under this hypothesis (to be discussed in the examples), applying the symplectic change ([3]) to ^ provides a generating 
function, that we denote S s (q), for the global manifold W s in the coordinates (q,p): 

P = X7S s (q), qeUnx^iU)- (6) 



It is easy to check, from the expression of the change ([3|), that S s = S s o \ + const. In this way, both manifolds W u 
and W s are expressed, in their common neighborhood 14 ("1 U, in terms of the same coordinates (q,p). Hence, we can 
study the transversality of their intersection along 7 by comparing the generating functions S u (q) and S s (q). 

More precisely, the transversality of the invariant manifolds will be studied by comparing, inside the 3-dimensional 
energy level containing them, the slope of the coordinate P2 with respect q2, which is a transverse direction to 7. Since 

dS u dS s 

in ( 4 ) and ( 6 ) we have p 2 = ^ — an d Pi = -7; — respectively, we define the functions 

dq 2 oq 2 



Such functions are defined, respectively, in intervals UT\{q 2 = 0} and I :=UC\x 1 (^)n{<72 = 0}. Actually, a relation 
between T s (qi) and the function 

n ?1 ):=^,0) 

can be explicitly given. In terms of such functions, our main result can be stated as follows. 



Theorem 1 The functions T n (qi) and T s (q\) are solutions of two Riccati equations whose coefficients are given 
explicitly from the coefficients up to order 2 in the expansion with respect to q 2 of B(q) and V(q). Assuming that they 
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can be extended up to the common point referred to in hypothesis (H4), namely q\ or q\ — Xoi.1i) respectively, then 
a necessary and sufficient condition for the transversality of the invariant manifolds W u and W s along 7 is that the 
following inequality is fulfilled: 

T a (qt)*f'(qt), (8) 

where the function T s (qi) can be expressed, through the change Q), in terms ofT s (qi), B(qi,0) and V(qi,0). The 
transversality is kept along the whole orbit 7. 



This theorem is deduced in Section 2.4 where we provide explicitly the Riccati equation for the function T n (qi) 
(see Theorem [6]), and show that it has a unique solution under a suitable initial condition. To obtain this Riccati 
equation, the key point is to consider the generating function S u (q) as a solution of the Hamilton-,] acobi equation: 

H(q,VS"(q)) = 0, 

and to use the expansion in qi of this equation. As far as T n (qi) is bounded, the unstable manifold W u admits a 
generating function S u (q) as in Q. Of course, similar considerations can be formulated for the function T s (q~i), leading 
to an analogous Riccati equation. Nevertheless, many examples satisfy a reversibility condition and the solution of the 



Riccati equation for the stable case can be deduced directly from the unstable one (see Section 2.5). The transversality 
condition ([8| can be checked from the Riccati equation, even in some cases where it cannot be solved explicitly, through 
a qualitative study of its phase portrait, as we show in some examples. 

Organization of the paper. The results of Section [2] are valid for any 2-dimensional configuration manifold Q, 
whereas in Sections [3] and [4] we restrict ourselves to the case of a 2-torus: Q = T 2 . On the other hand, we deal with a 
non-perturbative situation in Sections [2] and [3j and with a perturbative situation in Section |4j 

We start in Section [2~T] by studying the invariant manifolds W u ' s in a neighborhood of the hyperbolic point O, 



showing that they can be expressed in terms of (local) generating functions S n ' s (q). In Section 2.2 we consider the 
extension of the generating function S n (q) along the homoclinic/heteroclinic orbit 7, and consider its Taylor expansion 
in a transverse direction. We also discuss the transversality condition, from the comparison of the generating functions 
S u (q) and S s (q). In fact, we give the results for the unstable manifold of the point O, whereas the results for the 
stable manifold of O would be analogous. Expanding in q 2 the generating function S u (q), we see that the coefficients 
of orders and 1 are given by the orbit 7, and the coefficient of order 2 is the function T u (qi), which allows us to 
give a necessary and sufficient condition for the transversality of the invariant manifolds along 7. The fact that the 
outgoing part of 7 is contained in qi — imposes several restrictions on the Hamiltonian and the generating functions, 
which are studied in Section 2.3 Next, we deduce in Section 2.4 the Riccati equation (Theorem [6]) by using that 



the generating functions are solutions of the Hamilton- Jacobi equation, and we also establish the appropiate initial 
condition for the Riccati equation. As pointed out above, one should consider two Riccati equations, for both the 
unstable and the stable manifolds. Nevertheless, we show in Section [275] that, if we consider a reversible Hamiltonian, 
the solution of the Riccati equation for the stable manifold can be deduced from the unstable one, and in this way it 
is enough to consider one Riccati equation. As an example, in Section [2^6] we revisit a Neumann's problem considered 
in |Dev78j . consisting of an integrable system on Q = S 2 with transversality of the invariant manifolds. 

In Section |3.1[ we reformulate the statements and results of Section [2] for the case of the configuration manifold 
Q = T 2 , where we can take advantage of the periodicity in q±, and 7 is a homoclinic orbit or loop. This case is 
interesting in view of its relation to double resonances in nearly-integrable Hamiltonians with more than 2 degrees of 



freedom (see Section 1.2 1. As an application, we consider in Section 3.2 the example of two identical pendula connected 
by an interacting potential, generalizing the results obtained in GS95 for the case of a linear interaction. Using bounds 
of the solution of the Riccati equation for this case, we provide a sufficient condition for the transversality. 



The last part is devoted to the case of a perturbed Hamiltonian H e = H + eH*. In Section 4.1 we assume that, 
for £ = 0, the unperturbed Hamiltonian H has a loop 7 contained in qi — as in the previous sections. We consider 
two cases: (a) the unperturbed invariant manifolds intersect transversely along 7; and (b) the unperturbed invariant 
manifolds coincide, becoming a 2-dimensional separatrix containing 7 as an orbit belonging to a 1-parameter family 
of loops 7 S , s € K. In Theorem |16| we provide sufficient conditions ensuring, for e 7^ small enough, the existence of 
a perturbed loop 7 £ and the perturbed invariant manifolds intersect transversely along it, although in the case (b) we 
have to impose an additional condition on the perturbation iJ* (we stress that the existence of a loop is assumed in 



the unperturbed Hamiltonian H, but not in the perturbed one H £ ). Next, we show in Theorem 17 that the additional 
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condition for the case (b) can be expressed in terms of a Mel'nikov potential, defined in (81 -85) as the integral of 
along the unperturbed loops j s 



which allows one to check the transversality in concrete cases. Finally, in Section [4.2| 
we apply these results to the an example consisting of two pendula with different Lyapunov exponents, plus a small 
interacting potential of order e. 



1.2 Motivation 



The main interest for the 2-d.o.f. model considered in this paper lies in its close relation to resonant normal forms. For 
a given nearly-integrable Hamiltonian H(ip,I) — h(I) + sf(ip, I), with N > 2 degrees of freedom, in action-angle vari- 
ables, the mechanism described in |Arn64j to detect instability (Arnol'd diffusion) is based on the connections between 
invariant manifolds of (N — l)-dimcnsional hyperbolic invariant tori, associated to simple resonances. Nevertheless, 
along the simple resonances one also finds double resonances, which should be taken into account. 

Let us give a brief description of resonant normal forms in this context (for details, see for instance |BG86| . [LMS03 , 
ch. 2] and also |DG01j ; the ideas were initially developed in [Nck77 ). To study the behavior of the trajectories of 
H in the region close to a resonance of multiplicity n, with 1 < n < N (associated to a module of resonances 
jVf C Z N ), one carries out some steps of normalizing transformation, in order to minimize the nonresonant terms of 
the Fourier expansion in the angular variables ip. In this way, one obtains a symplectic transformation $ leading to new 
variables (ip, J), in which the Hamiltonian becomes a resonant normal form plus a small remainder: % o $ = T + R, 
where T only depends on resonant combinations of angles. By means of a linear change, we can assume that T 
only depends on . . . , j> n ) = q. Writing tji = {q,ip) <G T" x T N ~ n and J = (p,J) e R" x R N ~ n , we have 
r = T(q,p, J), R = R(q,ip,p, J), and we can study the (truncated) normal form T as a first approximation for the 
whole Hamiltonian T-L. 

If we neglect the remainder, for the normal form T(q,p, J) we have J = 0. Then, taking J as a parameter we can 
consider an n-d.o.f. reduced normal form, Fj(g,p) = T(q,p,J), and the behavior in the coordinates (tp,J) becomes 

just a set of N — n rotors: ip = J = 0. It is very important to understand the behavior of the reduced normal 
form in the coordinates (q,p) because this gives a first approximation to the original N-d.o.i. Hamiltonian T-L. 

Under certain conditions, the reduced normal form T°j(q,p) has equilibrium points, which provide a first approxima- 
tion for resonant (N — n) -dimensional invariant tori of the whole Hamiltonian %. In the same way, the n-dimensional 
invariant manifolds of hyperbolic equilibrium points of the reduced normal form, provide approximations for the 
./V-dimensional invariant manifolds associated to hyperbolic resonant tori, also called hyperbolic KAM tori (whose 
splitting seems to be closely related to Arnol'd diffusion). 

We point out that the "reduced" invariant manifolds can easily be studied in the case of a simple resonance 
(n = I), because in this case the I-d.o.f. reduced normal form is integrable and the invariant manifolds become 
generically homoclinic connections. For the whole Hamiltonian H, the existence of homoclinic intersections between 
the Af-dimensional manifolds of hyperbolic KAM tori was shown in |Eli94j , [DGOOj . Besides, the Poincare-Mel'nikov 
method can be used to measure the splitting of the invariant manifolds in some restrictive models (see [LMS03 , 
|DG04| for N = 3). As another related situation, the case of a loop asymptotic to a center-center-saddle equilibrium 
was considered in KLDG05 , proving the existence of homoclinic intersections and their transversality for hyperbolic 
KAM tori, contained in the center manifold and close to the equilibrium point. 

But for a multiple resonance (n > 2), in general the n-d.o.f. reduced normal form is non-integrable and hence 
the behavior of its invariant manifolds cannot be fully understood, although it is possible to give some partial results 
which can be useful for the whole Hamiltonian. In this context, it is proved in [LMS03, §1.10.2] that, if there exists a 
homoclinic orbit 7 for the reduced normal form such that the invariant manifolds W u,s intersect transversely along 7, 
then one can establish the existence of intersection for the invariant manifolds of the whole Hamiltonian %, together 
with a lower bound (N — n + 1) for the number of homoclinic orbits. A similar result is proved in |RT06] . for a 
model for the behavior near multiple resonances. Concerning a double resonance (n = 2), the case of a Hamiltonian 
with a loop asymptotic to an equilibrium point with 2 saddles and N — 2 centers has recently been considered in 
|DGKP10) . showing under some restrictions the effective existence of (transverse) homoclinic intersections associated 
to hyperbolic KAM tori on the center manifold. On the other hand, the dynamics near a double resonance has been 
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studied in |Hal95j and |Hal97) . but assuming that one of the resonances is strong and the other one is weak (in this 
situation, the reduced system is close to integrable). 



Coming back to the classical model considered in Section [Llj the Hamiltonian H defined in ([I]) can be seen as a 

particular model for the reduced normal form, where our assumption that B(q) is positive definite can be related with 
a quasiconvexity condition on the integrable part h(I), and the potential V(q) comes from the perturbation f(ip,I) 
via the resonant normal form. For such a model, we are assuming that we are able to describe a homoclinic orbit 
7, and we are going to provide a condition which allows to check the transversality of the invariant manifolds W u ' s 
along 7 in concrete examples. 

We point out that the existence of homoclinic orbits for a Hamiltonian of type ([!]) can be established using 
variational methods. It is shown in |Bol78j (see also |BR98a| . |BR98b| ) that, if the configuration manifold Q is 
compact (such as Q = T 2 ) and the potential V(q) has a unique maximum point, which is nondegenerate, then there 
exists a homoclinic orbit to this point. 



2 The generating functions as solutions of a Riccati equation 

2.1 The local generating functions around a hyperbolic equilibrium point 

We show in this section that, in a neighbourhood IA of a hyperbolic equilibrium point O, the local invariant manifolds 
W u ' s can be written in terms of generating functions: p = WS u ' s (q). Although most of the results of this section 
are standard, their proof is included here for the sake of completeness and notational convenience. At the end of the 
section, we also establish some local properties of the functions S u,s (q). 

To carry out this local study, we use the quadratic part of the Hamiltonian function (fTl) : 



where A = — T) 2 V(0, 0) is a symmetric matrix, and B = Bg(0) = -8(0, 0) is a positive definite symmetric matrix. First 
we show that, under hypothesis (HI), the point O at the origin of the coordinates (q,p) is a hyperbolic equilibrium. 

Proposition 2 If the potential V(q) of the Hamiltonian has a nondegenerate maximum at the origin q — (0, 0) of 
the configuration space, then (q,p) = (0,0,0,0) is an equilibrium point of hyperbolic type of the Hamiltonian. 

The following lemma will be used in the proof of this proposition. 

Lemma 3 Let G and Q be two real symmetric matrices with G positive definite. Then, the product GQ has all its 
eigenvalues real and positive if and only if Q is positive definite. 

Proof of Lemma\^ Indeed, let us consider the Cholesky decomposition G = LL T . It turns out that GQ and L T QL 
are similar matrices, since GQ = LL T Q ~ L- l LL T QL = L T QL. Therefore, the lemma follows at once from the 
application of Sylvester's law of inertia. □ 

Proof of Proposition^ If the potential V(q) has a nondegenerate maximum at q — (0,0), then the matrix A = 
— D 2 y(0, 0) is positive definite. The Lyapunov exponents of the equilibrium point O are the eigenvalues of the 
differential matrix of the Hamiltonian vector field: 




H{q,p) = \ (Bp,p) - | (Aq,q) + 0(\(q,p)\ 3 ), 



(9) 




To show that these eigenvalues are all real, we consider the square, 




G 



where we have (BA) T = AB, since A and B are both symmetric matrices. Then, by Lemma [3] all the eigenvalues of 
BA are real and positive, for A is positive definite. We also see that that the product matrix BA diagonalizes. Indeed, 
if we consider the Cholesky decomposition B = LL T , we see that L~ X BAL = L T AL, which is a symmetric matrix 
and diagonalizes, and so does BA. 

Let A- 2 , A| be the eigenvalues of BA, and V\,v 2 the corresponding eigenvectors. Since AB is the transpose of BA, 
it also diagonalizes with the same eigenvalues, whereas one sees immediately that w% — B~ 1 v\, w 2 = B~ 1 v 2 constitute 
a basis of eigenvectors: 

ABw k = Av k = B~ 1 BAv k = B^ 1 ■ \%v k = \£w k , k = 1, 2. 
Finally, the matrix DXh(0) has ±Ai,±A2 as eigenvalues with the four associated eigenvectors 



vi \ ± _ / v 2 

±AiWi ) ' 2 I ±A 2 w 2 



(10) 



respectively, as can be easily checked: 



±AfcUfe \ . ± 
X 2 k w k )= ±X ^- 

Thus, the origin (q,p) = (0,0) is an equilibrium point of hyperbolic type with Lyapunov exponents ±Ai, ±A 2 (real 
and nonvanishing) , and the proof of Proposition [2] is complete. □ 

Existence of generating functions for the (local) invariant manifolds. We have shown in the previous 
proof that the matrix BA diagonalizes, with real and positive eigenvalues. Hence, we can write 

BA = MA 2 lvr 1 , AB = NA 2 N~ 1 , 

where we define 

A:=^ Al ^ ^J, M:={vi v 2 ) , N := ( Wl w 2 )=B- 1 M. (11) 
We can assume, without loss of generality, that both Ai and X 2 are positive. 



The eigenvectors w^,w 2 of DXjjiO) introduced in (10) give the linear approximation of the (local) invariant 
manifolds W u,s . To be more precise, the eigenvectors , w 2 (with positive eigenvalues) give the unstable manifold 
W u , and the eigenvectors w^, w 2 (with negative eigenvalues) give the stable manifold W u . Hence, up to first order, 
each local invariant manifold can be parameterized as 

p = NA( + 0(\(\ 2 ), { p = -NA( + 0(\(\ 2 ), 

with q = (qi,q 2 ), p = (pi,P2), and parameters C = (Ci , C2) ■ On the other hand, since detM ^ 0, then by the implicit 
function theorem: 

C = M- 1 q + 0(\q\ 2 ), 



locally, in a neighborhood of q = (0,0). Therefore, substitution into the second equations of (12 1 yields, 

p = g"(q) = NAM- 1 q + 0(\q\ 2 ), p = g s (q) = -NAM-\ + 0(\q\ 2 ). (13) 

Thus, the manifolds W u,s can be expressed locally as the graphs. Moreover, due to the fact that W u ' s are Lagrangian 
manifolds, the restriction of the standard 1-form 6 — p\ dqi + p 2 dq 2 on W u ' s , 

Qlw^ = 9i' s dqi + g 2 s dq 2 , 



is a closed 1-form, and hence locally exact. Then, we have in ( 13 1 that the expressions g u,s = (g" ,s , g 2 ,s ) are the gradients 
of generating functions S u,s defined in a neighborhood of q = (0, 0) (and uniquely determined up to constants). Besides, 
it is not hard to see that the generating functions are as smooth (C r or analytic) as the initial Hamiltonian H. 
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Remark. The result given above is local, but suitable for our purposes since, as far as the invariant manifolds W u ' s can 
be expressed as graphs p = g n ' s (q), they admit generating functions of the form S n ' s (q) beyond a small neighborhood 
of the origin. Actually, it is well-known that the manifolds W u,s , which are asymptotic to an equilibrium point, are 
exact Lagrangian manifolds, that is #|yy u . a = dS' u,s with generating functions defined globally in the whole manifolds, 
gu,s . yyu,s — y ^ gee £ or j ns t ance |DR97j . |LMS03| ). Nevertheless, the global manifolds are not necessarily graphs 
P = P{o)i an d m such a case the generating functions have to be expressed in some other variables. 



Now, we define the following symmetric 2x2 matrices 

f u,s . = D 2 5 u,s( 0)0 ) 

The next lemma states some useful relations between the matrices A, B and £ u ' : 



(14) 



Lemma 4 For the matrices £ n ' s , we have: 

(a) £ U B£ U = A; 

(b) £ u = NAM^ 1 , positive definite; 

(c) £ s = —£ u , negative definite. 



Proof. We show that part (a) comes from the expansion of the Hamilton- Jacobi equation in q, taking the second-order 
terms. Indeed, expanding the gradient p — VS u ' s (q) at q = we have 



(15) 



VS u - s {q) = V i^S u s {0) + VS u ' s (0)q + - (D 2 S u ' s (0)g, q) + 0(\q\ 3 ) 
= VS u - s {0) + D 2 S u ' s (0)<? + OQq\ 2 ) = £ u ' s q + 0{\q\ 2 ) 



(recall that VS' u ' s (0) = for the origin of the coordinates (q,p) belongs to both W u,s , so no constant terms might 
appear in their expansions, compare also (13l). Next, substitution of ( 15 ) in ^ leads to the following second-order 
expansion of the Hamilton-Jacobi equation: 

H(q,VS»(q)) = H (q,£»q + 0(\q\ 2 )) 

= \ (B (£»q + 0(\q\ 2 )) , rq + 0(\q\ 2 )) \ {Aq, q) + 0(|g| 3 ) 

= ^(B£ u q,£ u q)-±(Aq,q)+O(\q\ 3 )=0, 
from which the desired equality of part (a) follows at once. 



Let us show (b). The expansions (13) and (15) give W U:S , locally, as the graph of a function; therefore we can 
identify 

£ U = NAM~ 1 , £ s = -iVAiLT 1 . (16) 



We see from the first of (16) that B£ u = BNkM~ l = MAM -1 (see the definition of matrix N in (11)), hence 



B£ u ~ A, which has positive eigenvalues, and B is positive definite. By Lemma|3] the matrix £ u is positive definite as 
well, proving part (b). Finally, we deduce (c) as an immediate consequence of the second of (16). □ 



Remark. Although the results of this section have been stated for the case of 2 degrees of freedom — so the matrices 
they refer to are 2x2 matrices — , the same results apply for n-degree-of-freedom Hamiltonians and consequently, for 
n x n matrices. 
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2.2 The generating functions around a homoclinic or heteroclinic orbit 



So far we have shown the local existence of the generating functions S u ' s (q) in a small neighborhood of the hyperbolic 
point O, as well as the tangent planes at O of the local unstable and stable manifolds W u,s . Now, our aim is to study 
whether such generating functions can be continued along the homoclinic/heteroclinic orbit 7, in order to study their 
transversality. Recall from hypothesis (H2) that, in the coordinates (q,p) of the neighborhood U, we are assuming 
that this orbit is contained in (72 = 0. In fact, we may consider such expansions for any Lagrangian manifold containing 
an orbit satisfying (H2), regardless of the fact that the orbit is asymptotic to hyperbolic points. 

Expansion of the generating functions in a transverse direction. To start, we consider the expansion in q 2 
of the potential V(q) and the matrix B{q) of the Hamiltonian ([lj: 

V(q) = V (qi) + V 1 (qi)q 2 -±Y(q 1 )qi + O(q2), (17) 
B(q) = B (q 1 )+ 1 ,B 2 (q 1 )qi + O(qi), (18) 



ith 



B ^ qi > = I h ( \ h f \ ' 3 = > 2 



(recall that we have Bi(qi) = according to hypothesis (H3)). Notice that the matrices introduced in ([9| can be 
expressed in terms of such expansions: 

.( -C(0) -V{{0) \ / 6 U0 (0) 6i 2 o(0) \ 

{ -V{(0) Y(Q) J ' 1 6x20(0) 6220(0) j ^ 



Here, we study the function S n (q) near the hyperbolic equilibrium point O, as well as its continuation along 7, for 
q 2 close to and (qi,q 2 ) € U, with the help of the Taylor expansion in the variable q% (notice that we are using qi to 
parameterize 7 inside the neighborhood U, and qi provides a transverse direction to it). For both generating functions 
S u ' s (q), associated to the invariant manifolds W u ' s of the point O, we consider the expansions 

S^ s (q) = S^ s ( qi ) + S^iqi) 92 + ^(qi) ? 2 2 + O(qi). (19) 

The term of order in these expansions is determined up to an additive constant that, to fix ideas, will be chosen in 
such a way that Sq' s (0) = 0. Notice that the matrices defined in (14) are 



(^• s )"(o) (sn'(o) 



(20) 



^ (Si'T(0) T u > s (0) J ' 

Analogously, we consider the generating functions for the invariant manifolds W u ' s of the point O, 

S^ s (q) = + S?' s (Qi) q 2 + |f ^(q,) qi + 0(q$), (21) 

also with Sq' (0) = (we show in the subsequent sections that, if a suitable symmetry occurs, the generating functions 
around O can be deduced from the ones obtained for O). 



Under hypothesis (H4) , there is a common neighborhood IA PiU where we can apply the change |3| , induced by 
q = x{q)- With this change, we can write (a piece of) the stable manifold W s of O in the coordinates (q,p), in order 
to compare it with with the unstable manifold W u of O. Applying this change to equation ([5]), we obtain the equation 
p = Dx(<7) T V5 s (x(g)) = V(S' S o x)(g), i.e. equation M, also with a generating function 



S s {q) = S s o X {q) + const (22) 

(in this case, the additive constant will not be taken equal to zero; see below). Now, we can consider an analogous 
expansion of the function S s (q), for q 2 close to and {qi,q 2 ) & U n X~ 1 {U), 

S s (q) = SS( ft ) + Sl(qi) 92 + ^T s ( qi ) qi + 0(qi), (23) 
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where the coefficients can be determined from the ones in (21 1, applying the change x- In particular, th e fu nction 



T s (qi) can be determined from 5q(?i), <Si(<7i) and T s (qi) (for an illustration, see the example of Section 2.6). Our 
aim is to compare the expansions of S u (q) and S s (q) in their common domain. 



It is important to stress that the coefficients of orders and 1 in the expansions ( 19 ) and ([23]) are determined by 
the orbit 7 itself. Indeed, if we consider the function S u (q), we see from Q and hypothesis (H2) that the orbit 7 is 
given, in the neighborhood U 1 by the equations 



7 



Pi 



dS u 
dqi 



P2 



dS u 



(gi,0) = S?(<zi) 



(24) 



(notice that at the hyperbolic point O we have (Sq)'(O) = S"(0) = 0). Since the same can be done with the function 
S s (q), we deduce that the coefficients of orders and 1 for both functions coincide. We can introduce a common 
notation for them: 

SoOzi) := S'o(Qi) = S s (qi), Si(<&) := = S!(«i), (25) 

for any qi £ I, the common interval for both generating functions. Notice that we can choose the additive constant 



in (22 ) in such a way no additive constant appears in the first equality of (25 1 (this is not very important because the 



constant does not take part in the gradient equations, but it is useful in order to fix ideas). 



The transversality condition. Next, with the help of the expansions introduced above we provide the condition 
for the transversality of the 2-dimensional invariant manifolds W u and W s along 7. This transversality must be 
considered inside the 3-dimensional energy level Af containing both invariant manifolds. Since Af is given by the 
equation 11 = 0, and on 7 we have = Qi 0, by the implicit function theorem we have that, near 7, the energy 
level Af can be parameterized as pi = g(qi, 92, P2) ■ In the coordinates (qi,q2,P2) 01 Af, from Q we see that the 
unstable manifold W u is given by the equation 

P2 = ^-(«i,8a) = SM + T"( qi ) q 2 + 0(qi). 
oqi 

This says that, in the coordinates (91,92,^2), the coefficient T u (qi) provides the slope of the manifold W u in the 
direction of q 2 , which is transverse to 7. Analogous considerations can be formulated for the stable manifold W s , 
whose slope is given by T s (qi). Then, a necessary and sufficient condition for the transversality of W u and W s along 7 
is that the two slopes are different for some fixed q\ £ I, as stated in pi). 



2.3 Restrictions on the functions defining the Hamiltonian 



We assumed in hypothesis (H2) that the orbit 7 satisfies q 2 = in the neighborhood U where the coordinates (q,p) 
arc defined. In this section, we show that this hypothesis implies several equalities, that will be used later, involving 



the coefficients Vb(<7i), Vi(gi), (Sq)'(<7i) an d Si(qi), appearing in the expansions (17) and (19) 



First of all, we parameterize the ourgoing part of the orbit 7 (inside the neighbourhood IA) as a trajectory. According 
to hypothesis (H2), we have 



31 = «?(*), 92 = 0, pi=p?(t), P2=P° 2 (t), t<h, 



(26) 



where 3i(t) > is an increasing function asymptotic to as t — > —00. Reparameterizing 7 as a function of the 
coordinate qi, we obtain the functions (Sq)' and S 1 " as in (24). Recall that, to simplify the notation, in (25) we have 
rewritten those functions as S' and S\ respectively. 



We also define the function 



dot gp(gi) 
^220(91) 



(27) 



Proposition 5 

(a) The inner dynamics along 7 in the neighborhood U is given by the differential equation q\ = [3(qi)S' (qi). 
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(b) The functions in (25) are given by 



-2V ( qi ) &i2o(gi) ql( v 

V PWl) 2 2o(9l) 

(c) TTie following equality is satisfied: S[(qi) fl(qi) S' (qi) = —Vi(qi). 

Proof. We consider the Hamiltonian equations ([2]), restricted to the orbit 7, as well as the fact that 7 is contained in 
the zero energy level of the Hamiltonian: 



qi = hw(qi)Pi + &i2o(<?i)P2, 
= &12o(<7l)Pl + b 2 2o{qi)p2, 
p\ = -±{B' ( qi )p,p)-Vt(qi), 
i>2 = -Vi(qi), 

H(q 1 ,0,p 1 ,p 2 ) = \ (bu (qi)p? + 2b 120 (qi)piP2 + b 22 o{qi)p 2 ) + V { qi ) = 0. 



(28) 
(29) 

(30) 
(31) 



According to (24), we can replace p\ — S' (q\) and P2 = Si(5i) 



definition of f3{qi) in (27), we obtain (a): 

<7i = ( ^110(91) 



As a direct consequence of (29), we obtain the second equality of (b). Replacing it in (28) and recalling the 

fri2o(<?i) 2 ' 



&220 (<?i) 



S' ( qi ) = (3( qi )S' ( qi ). 



Since qf (t) is increasing and /3(gi) > 0, we deduce that S' ( q i) > for <ji > 0. In the same way, we see that (31 1 can 
be written in the form 

i/3( gi )^( 9l ) 2 + K)(gi) = 0, 



which gives the first equality of (b). Finally, replacing p 2 — S[(qi) q± in (30), we obtain (c) 



□ 



Remark. In (b), we can write both Sq(c/i) and Si(qx) in terms of Vo(qi) and Bo(qi). Inserting this in (c), we obtain 
an equality allowing us to obtain an explicit expression for V\_(qi). In other words, the functions Vo(?i) and Vi(qi) 
in (17) cannot be independent, due to the existence of the orbit 7 that satisfies (J2 = 0. 



Remark. To obtain (30), we used that -Bi(gi) = in (18), according to hypothesis (H3). If this hypothesis is not 
assumed, after straightforward computations there appear more terms in the right hand side of (30) and, consequently, 
in the formula of item (c) of this proposition. 



2.4 The Riccati equation 



In this section, we show that the function T u ( q i), defined in Q or (19) from the generating function S u ( q ) of the 
unstable invariant manifold W u of the point O, is a solution of a first-order differential equation of Riccati type in 
the variable q%. This could be done in the same way for the function T s (qi) associated to the stable manifold W s 
of O, obtaining an analogous Riccati equation in q±. In this way, the solutions of such Riccati equations provide the 
slopes of the two invariant manifolds in a transverse direction to the orbit 7. We are assuming in our hypothesis (H4) 
that both solutions can be extended up to some cfc and g*, related by the change x in (|3|. This changes provides the 
relation between the expansions (21 ) and (23 1, allowing us to obtain the value of T s (qi) and compare it with T n (qi) in 
view of the transversality condition Thus, in principle both equations should be solved in order to decide whether 
the invariant manifolds are transverse along 7. 



Nevertheless, in many cases some kind of reversibility relations are fulfilled, with an involution relating the two 



invariant manifolds, as well as their generating functions (see Section 2.5), and it will be enough to find the solution 
of only one of the Riccati equations. Since the examples considered in this paper satisfy some reversibility, in this 
section we only deal with the function T u (qi) associated to the unstable manifold of O. 
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In order to formulate the Riccati equation for T u (qi), we define the functions 



S(qi) := &i2o(<7i)<S'i(<2 , i), 

a(qi) := Y(q x ) - bno{Qi)S[ (<?i) 2 

-\ [b ll2 {qi)S' Q { qi ) 2 + 2b 122 (q 1 )S' (q 1 )S 1 (q 1 ) + b 222 {q l )S l {q l f 



(32) 
(33) 



and recall that f3{qi) was defined in (127 



Theorem 6 The junction T u (qi) in the expansion {19) of the generating function of the unstable manifold W u ; is a 
solution of the Riccati equation: 



/3(qi)S' ( qi ) (T u )' + 26( qi ) T u + b 220 ( qi ) (T u ) 2 = a( Ql ) 



(34) 



Proof. We use that the generating function S u (q) is a solution of the Hamilton- Jacobi equation. Restricting the 
Hamiltonian ([I]) to the unstable manifold W u , we have an expansion 



H(q, VS u (q)) = H^ gi ) + HUqi)q 2 + W^q 2 + 0{ql) = 0, 



and we can replace 



V5 u (g) = 



qi 



where S^(q±) corresponds to the coefficient of third order in q 2 , in the expansion (19). Using also (17-18), we obtain 

K(qi) = \ [W<Zi)S (9i) 2 + 2b 120 (q 1 )S' (q 1 )S 1 (q 1 ) + 6 220 (gi)5i(<Zi) 2 ] 
+V (qi), 

Hl(qx) = b 110 (q 1 )S' (q 1 )S[(q 1 ) + b 120 ( qi ) lS' Q ( qi )T u ( Ql ) + SMS'M] 
+b 220 (qi)S 1 (q 1 )T u {q 1 ) + V 1 (q 1 ), 

H u 2 (qi) = b 11 o(q 1 )[S' (q 1 )(T^'(q 1 ) + S' 1 (q 1 f] 

+26i 2 o(gi) [±S' ( qi )S$( qi ) + S{(qi)T"( qi ) + ^faXT^'fa)] 
+6220(91) [Si(qi)S%( qi )+T»( qi ) 2 ] 

+ § [hl2(qi)S' {qx) 2 + 2b 122 (q 1 )S' (q 1 )S 1 (q 1 ) + b^q^S^q,) 2 ] 
-Y(qi). 

Since H(q, VS u (q)) = 0, it follows that Hj(qi) = 0, for j = 0, 1, 2. In particular, the expression of 'H 2 {qi) above yields, 
after some arrangements, 

[bno(qi)S' ( qi ) + b 120 (q 1 )S 1 (q 1 )} (T u )' + 2b 12() (q 1 )S , 1 (q 1 )T u + b 220 ( qi )(T u ) 2 
= Y(q 1 )-b 110 (q 1 )S[(q 1 ) 2 

- I [bn 2 (qi)S' ( qi ) 2 + 2b 122 (q 1 )S' (q 1 )S 1 (q 1 ) + b 222 (q 1 )S 1 ( 9l ) 2 ] , (35) 

where the terms multiplied by S^(qi) vanish, due to the equality ( [29] ) in the proof of Proposition [5] Now, taking into 
account (28) and Proposition [5] a), one has 

&110 (qi)S' (qi) + bi 20 (qi)Si(qi) = /3(qi)S' (qi) 



which, together with the definitions in ([27] of S(qi) and a(qi), gives rise to the Riccati equation (34) after substitution 
in p5l □ 



Remark. Recall that hypothesis (H3) is assumed throughout the computations. Otherwise, one reaches the same 



form (34) of the Riccati equation, but its coefficients S(qi) and a(q\) hold additional terms coming from the linear 



part -Bi(gi) of the expansion (18) in 172, of the matrix B(q) 
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The proof of Theorem [6] almost completes the proof of Theorem [7] To finish it, it is eno ugh to formulate the 
analogous Riccati equation for T s (q\), and take into account the considerations of Section 2.2 on the transversality 



condition ([8|), and on the relation between the functions T s (q~i) and T u (qi) (the latter one taking part in (|8j)). 
Concerning this relation, recall that the function T s (qi) can be determined from Sg(gi), S'f(gi) and T s (q~i). Now, 
according to the first remark after Proposition [5] (which applies also to the stable manifold), we see that T s (qi) can 
be determined from Vo(qi), £?o(<7i) an d T s (qi), which is the last assertion of Theorem [l] 



Coming again to the Riccati equation (34), it is clear that it has a singularity at q\ = since Sq(0) — 0. Thus, in 



principle the existence and uniqueness of solution might not take place. However, we are going to establish the right 
initial condition for the equation, and show that the solution is unique. 



Lemma 7 The function T n (qi) solving (34) satisfies the initial condition 

-5(0) + VA a(0) 



T u (0) = 



6220 (0) 



6(0) 



5(0) 2 + b 220 (0)a(0). 



(36) 



Proof. Notice that a solution of ( 34 ) defined at q\ — satisfies the equality 

28(0) T u (0) + 6220(0) T u (0) 2 = a(0), 

and hence the initial condition is almost determined by the differential equation: T u (0) = (—(5(0) ± vA)/&22o(0), 
where only the sign '±' has to be determined. 



Taking into account that <Sq(0) = S\(0) — 0, we deduce from Proposition |5jb) that 

Si(0) = -^@5?(0). 



We know from Lemma Qb) that the matrix £ u in (20) is positive definite. Then, we have S"(0) > and 

< S£(0)T u (0) - S[(0) 2 = Sft(0) 



T u (0) + ^og) 5 , (0) 



T u (0) 



6220(0) 
6(0) 



6220(0) 



and putting these formulas together, it turns out that we have to choose the positive sign in (36). 



□ 



Remark. Using similar arguments, it is also easy to show that A > 0. However, this is not necessary because the 



results of Section 2.1 ensure the existence of the generating function S u (q) around the origin and, subsequently, a real 
value for T u (0). 



Uniqueness of solution. Despite the singularity of (34) at q\ = 0, we show in the next proposition that there 



exists a unique solution satisfying the initial condition (36), which gives the function T u (q\) associated to the unstable 
manifold W n . 

Using the change of variable qi = q1(t), provided by the orbit 7 in the neighborhood we obtain another useful 
expression for the Riccati equation ( ]34[ ) , in terms of the time variable t. F or a ny given function f(qi), we write 
f(t) — f(qi(t)). Then, we see from Proposition [5^ a) that the Riccati equation (34) becomes 



f +28(t)T u + b 220 (t)(f u ) 2 =a(t). 



Proposition 8 The Riccati equation {34) has a unique solution that satisfies the initial condition (36) 
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Proof. Clearly, the unstable manifold W u exists, and its generating function provides a solution T u of the Riccati 
equation (34 1 with the initial condition (36 1. Let us show the uniqueness of solution. If T(qi) is another solution, the 
difference U(qi) — T(qi) — T u (qi) is a solution of the associated Bernoulli equation, with U(0) — 0. In terms of t, the 
Bernoulli equation for U(t) becomes 



U + ^(t)U + b 220 (t) u 2 = o, 



lim U(t) = 0, 



(37) 



where we define 



#) = 2(6(t) + b 220 (t)f u (t)). 



Denoting ipo = lim ip(t), we deduce from (36) that 

t— >— oo 



V^o - 2(5(0) + &22o(0)T u (0)) = 2VA > 



We prove the uniqueness of the solution U(t) = for (371 in a very simple way. The idea is that the linearization 
of (371 tends, as t — > — oo, to the equation U + ipoU — 0, with ijj > (i.e. the origin is unstable as t — > — oo). In fact, 
similar results for higher dimensions were established in |CL55I ch. 13]. Denoting u(t) = |[7(t) 2 > 0, we have: 

u(t) = U(t)U(t) = [-Vo ~ m) - <Ao) - b 2 2o(t)U(t)} ■ 2u(i). 
There exists t\ < such that \4){t) — r/'ol < ^0/2 and |&22o(t)^(i) < "0o/2 for any t < ii. Then, 

u(*) < [-V'o + V^o/2 + VV2] • 2u(t) = 0, t < ti, 
and we see that it is not possible to have lim u(t) = 0, unless u(t) = 0. □ 



Remark. One could expect that each solution of the Riccati equation (34 1 generates a 2-dimensional manifold 
consisting of a 1-parametric family of trajectories, containing 7. Only one of such manifolds, namely the one associated 
to the initial condition (36), is the unstable manifold W u . 



To end this section, we point out that one could use the variational equations around the orbit 7 as an alternative 
method in order to describe the invariant manifolds W u ' s . Such a method is followed (for a more particular case) 
in [GS95j. [RT06 . Since the variational equations are equivalent to a second-order linear differential equation, the 
well-known relation between such linear equations and Riccati equations via a change of variables provides a relation 
between the approach using the variational equations and our approach using the Hamilton-Jacobi equation, which 
leads to a Riccati equation. 



The advantatge of using the Riccati equation is that a qualitative analysis of its phase portrait can be carried out 
with the help of simple methods of dynamical systems. Such a qualitative approach is useful when the solutions of 



the variational equations cannot be obtained explicitly, or they have complicated expressions. In Sections 2.6 and 3.2 



we illustrate with some examples the use of the Riccati equation. First, in the example of Section |2.6| we show that 
this method is simpler than solving explicitly the corresponding second-order linear equation. In the example of 
Section [3~2} the linear equation can be solved explicitly in some particular cases, but in general it is not intcgrablc. 



2.5 Reversibility relations 



In this section, we assume that the Hamiltonian equations ([2]) satisfy a reversibility condition, which relates the two 
invariant manifolds W u,s of the hyperbolic point O. We are going to show that this reversibility implies a relation 
between the generating functions S u ' s (q) introduced in ( |19[ ), and hence between the associated slopes T u ' s (gi). In fact, 
we should consider the invariant manifolds W u and W s of the points O and O respectively, in order to study whether 
they are transverse along a co mmo n p iece of the orbit 7. Nevertheless, if the Hamiltonian satisfies some symmetry or 
some periodicity (see Sections 2.6 and 3.1| then it is not hard to relate the manifolds W s and W s . Consequently, we 



formulate the results of this section for the invariant manifolds W u,s , in the neighborhood U of the first point O. 
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It will be enough, for the examples considered, to consider the following type of reversibility. We say that H is 
IZ-reversible if the Hamiltonian equations ^ are reversible with respect to the linear involution 

TZ:(q,p)^(Rq,-Rp), (38) 

with a given matrix 

R = ^ n n,r 2 = ±l. (39) 

Lemma 9 The Hamiltonian H in |7p is IZ-reversible if the functions B(q) and V(q) satisfy the following identities: 

B(Rq) = RB(q) R, V(Rq) = V(q). 

Proof. It is well-known that the reversibility with respect to 7Z is equivalent to the identity Xh olZ = —TZXh- Using 
that 1Z T J7Z — — J (i.e. 1Z is 'antisymplectic'), the previous identity becomes V(H o TV) = VH, which can be written 
as H o 1Z = H + const. But the point O (the origin of the coordinates x = (q,p)) is a fixed point for 1Z, hence the 
functions H oTZ and H must coincide at this point and the constant vanishes. 

Now, we have the equality H(Rq, —Rp) = H(q,p). To finish the proof, it is enough to write H in terms of the 
functions B(q) and V(q) as in 0. □ 

Remark. This condition for the 7\L-reversibility of the Hamiltonian |l]) is always satisfied if we choose r% = r 2 = 1 



in (39), which can be called the trivial reversibility. This can be enough in some examples, but in other cases we will 



be interested in the other types of reversibility. 



Clearly, the reversibility gives a relation between the invariant manifolds of the hyperbolic point O: we have 
W s = TZ W u . In the next proposition, we deduce from this fact a relation between the generating functions associated 
to the invariant manifolds. 



Proposition 10 If the Hamiltonian H in is TZ-reversible, then the generating functions defined in (19) satisfy the 
identity 

S s (q) - -S u (Rq) 

and, for the coefficients of their expansions in qi, 



So(9i) = -SoCrift), S\{q x ) = -r 2 S?(r igi ), T s ( gi ) = -T u (r 1<Zl ). 



Proof. Recall that the manifold W u is given by the equation p — S7S u (q). Applying the reversibility (38), we 
obtain for W s the equation — Rp = \7S u (Rq), which must coincide with p = \7S s (q). Then, we have the equality 
VS S = —R VS U oR = -V(S U o R), which implies that S s = ~S U oR + const. Since in fHfy we set Sft(O) = S*g(0) = 0, 
it turns out that the constant vanishes. 



Expanding in q 2 and taking into account the form of the matrix R in ( 39 ) , we obtain as a simple consequence the 
equalities involving the functions Sg ,s (qi), Si' s (qi) and T u ' s (qi). □ 



2.6 Example: an integrable system on Q = S> 2 

We consider in this section the classical Neumann problem on the two 2-sphere S 2 (for an account see |Mos80j . and 
the same problem is tackled in |Dev78] . but working in the projective space P 2 ). This is an example of integrable 
Hamiltonian system with hyperbolic equilibrium points whose invariant manifolds intersect transversely along hete- 
roclinic orbits (or homoclinic orbits if one works in P 2 ). We are going to obtain this result of transversality applying 
our results on the generating functions of the invariant manifolds. 
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For 2 d.o.f., consider a particle moving on a 2-sphere S 2 in IR 3 under the action of the following potential on S 2 , 

V(x) = -±{Ax,Ax)=-~\Ax\\ (40) 

where A is a diagonal matrix, A — diag[Ai, A2, A3]. Further, it is assumed that A3 = 0, and that < Ai < A2. As 
usual, TS 2 denotes the tangent bundle, 



T§ 2 



Ux,x) £ I 3 x I 3 : \x\ 2 = 1, and (x,x) = o| 



Taking into account that VV(x) = —A 1 Ax = —A 2 x, the movement of particle on the sphere is described by the 
second-order equation in K 3 , 

x+ \x\ 2 x = A 2 x~ \Ax\ 2 x, (41) 

whose orbits lie on S 2 . It is not hard to check that this system has 2 hyperbolic equilibrium points at (0, 0, —1) and 
(0, 0, 1), with 4 heteroclinic orbits connecting them. In this section, our hcteroclinic orbit 7 will be the one that goes 
along the semi-circle X2 = 0, X\ > 0. 



The Hamiltonian. To tackle the equation (41 ), we introduce the two (local) charts: tp : R 2 — > U := § 2 \{(0, 0, 1)}, 
given by 

/ / \ ( 4gi 4 g2 q? + qj -4 \, 

9 = (91,92)^ v{qi,q2) = \t~, — 2~i — 2>r~i — 2~ — 2>r~i — 2~, — 2 > (42) 

\4 + q{+q£ b + q-f+qi 4 + q{ + q 2 2 ) 

and (p : R 2 — >U := § 2 \ {(0,0, -1)}, given by 

4<7i 4q 2 q? + q 2 2 



9 = (<?l»£fc) i-> ${q\Ai) = 



4 + q 2 + q 2 ' 4 + q 2 + q 2 ' 4 + q 2 + g 2 



so (U,(p) and (U,<p) are the usual stereographic coordinates in § 2 . These coordinates induce the corresponding two 
charts (U, ip*) and (U, (p*) in the cotangent bundle T*§ 2 , 

(q,p) £ I 2 x I 2 4 = (^( g ),D(^(g)- T p) eU ~ U x R 2 ; 

(q,p) e R 2 x M 2 h> <p*{q,p) = {^(q),V>(p{q)- T p) eU = UxR 2 . 

In particular, for the origin we have that (q,p) = (0,0) =: Oe R 2 x R 2 maps to 7> := p*(0,0) = (0,0, -1,0,0,0) in 
the chart (U,cp*) and (q,p) = (0,0) =: £ R 2 x R 2 maps to V := <p*(0,0) = (0,0,1,0,0,0), in the chart {U,(p*). 

On the other hand, we can give explicitly the change of coordinates or overlapping map between the two charts 
(U,tp) and (U,<p). This map X ■ S — > S, with £ := ^(i/ D t7) = M 2 \{(0,0)} and £ := ^(E/nE/) = R 2 \ {(0, 0)}, 
writes 

(9i,92) i-> (91,92) = x(9i,92) := o ^(91,92) = ( 2 ^ X 2 , tj-j] • (43) 

V 9i + 92 9i + 92 / 

Now, if we denote q = (91, 92) € S, q — (91,92) eS;p= (pi,p 2 ) € p = (pi, P2) € T|S, then, the extension of x 

to the cotangent bundles, x*, can be written as 

(q,p) e T*£ ^ ($,p) - X *(9,P) - (X(9),D X (9)" T P) € T*Z. (44) 



The system (41) is Lagrangian and, when expressed in the (intrinsic) coordinates (42), its Lagrangian function is 



L(q,q) = -(q,G(q)q)-V(q), 



(45) 



where G(q) = T)(p(q) T Dip(q) is the metric tensor of S 2 in the coordinates chosen, and V(q) = V o <p(q) corresponds to 
the restriction of the potential (40) on the 2-sphere. Explicitly, 



G(q) 



16 

(4 + 9! 2 + 9 2 2 f 



1 
1 



V(q) 



8(A 1 y + A 2 g2 2 ; 

(4 + 9! 2 + 92 2 ) 2 



(46) 
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Furthermore, the system (41) can be brought into Hamiltonian form, taking 



P = ipuPi) 



dL 



dq 



-(9,9) = G{q)q, 



as the actions conjugated to the coordinates q = (qi,q 2 ). Therefore, the corresponding Hamiltonian function writes 



where 



H(q,p) = ~{B(q)p,p)+V{q) 



B(q) = G(q)- 1 



(4 + ql + qt)(vl + gg) 8(A 2 g 2 + A 2 2 g|) 

32 (4 + 9l 2 + g|) 2 



(4 + g 2 +g 2 2 )' 
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We see from the quadratic part of (47) 



H 2 (q,p) 



1 
1 



p'i+p'i yfqf + xiq? 



(47) 



(48) 



that the origin (q,p) — (0,0), corresponding to the point V = (0,0,-1,0,0,0) in the chart (U,(p*) of T*§ , is a 



hyperbolic equilibrium point of the Hamiltonian flow associated to (41), with Lyapunov exponents ±Ai, ±A 2 . The 
same applies to V = (0,0, 1,0,0,0), which corresponds to the point (q,p) — (0,0) in the chart (U,(f*) of T*S 2 (see 
the remark below). We denote W u ' s the (local) unstable and stable manifolds of the point V and W u,s the (local) 
unstable and stable manifolds of the point V . 



Remark. The Lagrangian (45) takes the same form in either local coordinate system, (U,ip) or (U,ip). Thus, the 



associated Hamiltonian writes the same irrespectively of which coordinates (q,p) or (q,p) are taken in the phase space 
T*S 2 . Particularly, this implies that: (i) the point V, represented by (q,p) = (0,0) in the chart (14, tp*) is also a 
hyperbolic equilibrium point, and (ii) the (local) unstable and stable manifolds >V U ' S of the point V in the coordinates 
(q,p), and W u,s of the point V in the coordinates (q,p) are given, in terms of generating functions, by p = \7S n - s (q) 
and p — WS u ' s (q) respectively (i.e. we have S u,s = S u ' s as functions). 



Next, we expand the components of the matrix (48) in the Hamiltonian (47) in powers of q 2 , 

hj(q) = kjo(qi) + hji{q\)q2 + \b t0 2{qi)qi H — , i,j = 1,2. 

Notice that there are no terms of degree 1 in q 2 , i.e. biji(qi) = for i, j = 1, 2, which complies hypothesis (H3). For 
the remaining coefficients, we have: 



buo(qi 



&112(9l) 



(4 + gl 2 ) 
16 

4 + gf 
4 ' 



6l20(9l) = 6210(9l) = 0, & 22 ( 

6122(91) = 6212(91) = 0, 6222(91) 



(4 + 9ff 
16 ' 

4 + 9i 2 
4 



A similar expansion can be done for the restricted potential ( 46 ) 



with 



Vo(gi) 



V(q) = V (qi) + V 1 (q 1 )q 2 - \Y{ qi )q^ + 
8\?q? x » ,„ s 16 



(4+9i z r 



Vi(9i) = 0, F(9i) 



A 2 2Ajgf 
(4 + 9l 2 ) 2 V 2 4 + 9l 2 



The outgoing part of the heteroclinic orbit. Let us consider the local coordinates (U,<p). It can be checked 
out that the Hamiltonian (47) has as a solution a trajectory of the form |26)_[the outgoing part of 7) with q°(t) = 
(q1{t),q°(t)), p°(t) = (Pi(tfj%(t)) = B(q°(t))q°(t), where B(q) is the matrix 1 48}. Writting them down explicitly, 



qi = g°(t) = 2e^, ,J2 = 92 W = 0, 

16Aie Alt 



pi =mt) 



(4 + 4e 2A i*) 



r, P2=P° 2 (t)=0, 



(49) 
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which can be also parameterized by q\. 

7 : qi > 0, q 2 = 0, 



Pi 



16Aigi 
(4 + 9l 2 f 



P2 =0. 



(50) 



The corresponding trajectory on the phase space T*E> 2 connects the point V with the point V (as t goes from — oo 
to cxi, or qi goes from to oo). Furthermore, since q± increases along the orbit and q 2 = 0, the outgoing part of 7 
satisfies hypothesis (H2); then, according to (24 25) one has, for this example, 



pi = 5*0(91) = (<So)'(gi) 



16Ai<?i 
(4 + 9l 2 ) 2 



P2 = S 1 (q 1 )=S?(q 1 ) = 0, 



and the functions defined in ([27]) and ( |32fj33[ ) are found to be 

, 2 



%l)= (4 ^f ^ , 'M'/l > = n. n( 7l ) = 



1G 



16 



(4+ 9l 2 r 



A 9 - 



4A 1 y 

4 + ?1 2 



Remark. Note that, in agreement with Proposition [5](a), the outgoing part of 7 in (49) performs the inner dynamics 
qi = /3(gi)S , Q(gi) = \iq±, with the initial condition qi(0) = 2. 



Therefore, the Riccati equation (34) and its initial condition p6|, for the example at hand, are given by 

1G 



1 



Aigi(T u )'+-(4 + gi 2 ) z (T u ) 2 



(4+^r 



A 2 - 



4A 2 g 2 ' 
4 + ' 



T u (0) = A 2 . 



(51) 



Lemma 11 The solution of { 51), with the given initial condition, is defined for all q\ > 0. 



Proof. It is checked out immediately that 



I6A1 



(4+ 9l 2 r 



(52) 



is the solution of the Riccati equation (51) for A 2 = X%, with Tq'(0) = Ai. Let us denote by £(<7i) := T u (qi) — T^qi) 
and /1 := A 2 — Ai (/i > 0, since A2 > Ai); therefore, T u (q) is a solution of (51 1 iff £(gi) is a solution of 



A 1 g i r + 2A 1 e+-^(4 + g 2 ) 2 e 2 



1G 



16 (4 + 9! 2 ) 



(53) 



By Proposition |8j there exists e > 0, such that £(gi) is defined for < q\ < e. The idea is to extend this local solution. 
Consider < 77 < e such that £(77) > /a/2. From (53) it follows that £' > if qi > 77 and £ = 0, i.e. the direction field 
points upwards at £ = (see Figure 1(a) ); so £(<7i) > for q% > 77. On the other hand, 



£' = 



I6/A1 



( 2/J A 1+ /i 2 )--£- 

9i (4 + gf) 



qi ' I6A191 



< 



16 ^ (2fI W)-^ <4+ ' 2) ■= 



77(4 + ?7 2 )" 



91 



16Ai<?i 



for Qi > 77, with ci > 0. This means that the slope is bounded by a positive quantity and then, this local solution, 
£(<7i) with £(0) = (1, can be continued, and is defined for all qi > 0. □ 



Using the reversibility. One sees from Lemma [9] that the Hamiltonian ( 47 1 satisfies the reversibility defined in 
Section 2.5 with respect to the linear involution 1Z : (q,p) 1— > (-R<?, —Rp), with i? = diag[r!, r 2 ], and for any r 1; r 2 = ±1. 
In particular, if we take r± = r<i = 1, then Proposition [TO] yields 



2A o 2 

^(<Zi) = -5 u ( (Zl ) = -^^, <?f = -<?«(<?i) = 0, T s ( gi ) = -T u ( 9l ). 



(54) 
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Let us now consider the other local chart (U, (p*). The (local) stable manifolds W s of the origin O := (q,p) = (0, 0) 
(corresponding, in this coordinate system, to the point P — (0, 0, 1, 0, 0, 0) on the phase space T*S 2 ), can be put as a 
graph of a function through the same generating function 5 s , for theJ3amiltonian takes exactly the same expression 
in both local charts, as we pointed out before. Thus, the manifold W s is given by p = VS s (q). Therefore, we know 
from (22) that the change of coordinates (44) yields the identity S s = S s o \ + c, with c constant, and this gives W s 



in the coordinates (q,p) of the chart (U,(p). 



X = (Xi>Xz) S iven in < 43 >> one nas 



Expanding S s around the outgoing part of 7 given by (50), regarding the reversibility (54) and the components 



S s (q) = S* (q 1 ) + S s 1 (q 1 )q 2 + -f s (q 1 )qi + 



C-SZ( X l(q)) ~ S?(Xl(q))X2(q) ^ ( Xl (q)) X2 (q) 2 



c-ss 



4<?i 



q'i + ?! 



4qi 



9? + 92 J (9i 2 + 9 2 2 ) 



2\2 



and by comparison of coefficients, 



c-S% 



3i 



M 1 

Qi \9i 



( l2 



S s ( qi ) = c-SZ 



T s (qi) = ^(SW 
9l 



qi 



8A1 

"^T4' 

> (± 

Qi \Qi 



32A: _ 16 / £ 
fe 2 +4) 2 qf Ui 



(according to (25), the constant c can be fixed if, for example, we set Sq(2) = Sq(2) — Ai; this exacts c = 2Ai). Then, 

Ai 



the difference between T u and T s at q\ = 2 is 

T u (2)-f s (2) = 2 T U (2)- 
and therefore, 



T u (2) ^ 



Ai 



(55) 



1.5 



0.5 



111/77/ 

1117 7 7/ 



0.5 



1.5 



2.5 



(a) 




(b) 



Figure 1: (a) Direction field corresponding to the Riccati equation ,'53), where the solution with £(0) = /i := A2 — Ai, for A2 = 2, 
Ai = 1, is also plotted, (b) Three solutions of the Riccati equation \53\l for A2 =2; the lowermost one corresponds to Ai — 2 



and the other two, in ascending order, to Xi = 1.75 and Xi = 1, respectively. Note that the corresponding to Ai = A2 = 2, holds 
the point (2, Ai/4) = (2, 1/2). 
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is a necessary and sufficient condition for transversality between W u and W s along the heteroclinic orbit 7. This is 
stated in the proposition below. Notice that, in the chart (42 1, for q = (q*,0) = (2, 0) we obtain the point (1, 0, 0) € § 2 . 



Proposition 12 // in the system (41) defined on the 2-sphere S 2 
we assumej\2 > Ai > A3 = 
manifold W s of the point V = 



with A the diagonal matrix A — diag[Ai, A2, A3], 
then, for its associated Hamiltonian flow in the phase space T*S> 2 , the stable invariant 
(0, 0, 1, 0, 0, 0) and the unstable invariant manifold W u of the point V = (0, 0, —1, 0, 0, 0) 



intersect transversely along the orbit 7 (see (50)). 



Proof. By the considerations in the previous paragraph, it suffices to check the transversality condition (551, which is 

it is 



equivalent to £(2) ^ 0, where £(qi) := T u (qi) — Tq^i), with T^((?i) given by (52 1; but in the proof of Lemma 11 
stated that £ > for qi > 0. Thus, in particular, £(2) ^ is satisfied. This proves the proposition. 



□ 



Integrability of the Riccati equation. Alternatively, the Riccati equation (51 1 could have been solved explicitly 



in order to check the transversality condition (55). First, a standard change of type 

16Ai<?i 



y_ 
y 



(56) 



where T u — T n (qi), and y = y(qi) is the new unknown function, transforms (51 1 into the second-order linear equation 



3ft 2 



ft(4 + ft 2 ) 



4A 2 + (A 2 -4A 2 )ft 2 
A 1 V(4 + q 2 ) 



y = 0. 



(57) 



Remark. This is a particular case of the well-known fact that any Riccati equation z" + a(x)z + b(x)z 2 — c(x) 
transforms into a second-order linear equation, y" 



b' (x) \ y 
a(x) — I y' — b(x)c(x)y — 0, through the change z = 



b(x) 



b(x)y 



Now, the Kovacic 's algorithm can be applied to investigate the (Liouvillian) integrability of second-order differential 
equations with coefficients in the class C(x) of rational functions (for a complete description of the process and 
applications, we point the reader to |Kov86j . |DL92| . |AMWli] and references therein). This algorithm is implemented 
in some computer algebra systems such as Maple; and when applied to equation (57), it produces the following two 
fundamental solutions: 

Vi{qx) = (-4(A 2 + Ai) + (Ax - A 2 )q 2 ) q^ /X \ 
4(A 2 - Ai) + (Ai + \ 2 )q? — A2M1 

2/2(91) = r — t q x ; 



hence, the linear equation (57) is integrable and so it is the Riccati equation (51). Actually, transforming back 2/1(51) 

16Ai<?i 2/i(«i) 



using (56 1, one obtains a solution of \hl\ 
T"(qi 



(16 A| + 16 A1A2) (ft 2 + 4) + 32 \fq 2 



(4 + g 2 ) 2 yM ( g 2 +4 ) 2 ((Ai- A 2 )g 2 -4A 1 -4A 2 ) 



with T u (0) = A 2 ; and since A 2 > A x > 0, it is well defined for all q 1 > 0. Besides, the transversality condition (55) can 
also be checked straightforward from this explicit form of the solution of the Riccati equation (51 1 . Indeed, for qy = 2 
it gives 



T u (2) = i(A 2 + A 1 - 



A 2 



> 



Ai 



for A 2 > Ai > 0. In Figure |T(b) the function T u is plotted for different values of the parameters. 
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3 The case of a 2-torus: Q = T 2 



3.1 The generating functions in the periodic case 

In what comes, that the configuration space will be the 2-torus, i.e. Q = T 2 = (R/2irZ) 2 . In view of the form of 
the Hamiltonian (JlJ, this is equivalent to suppose that the quadratic form B(q) and the potential V(q) are both 
2ir-periodic functions with respect to q = ((71,(72)- In this way, the hyperbolic point O — (0,0,0,0) can be identified 
with O = (2-7T, 0,0,0), and therefore a biasymptotic orbit 7 connecting both ends is a homoclinic orbit or loop. 

Let us introduce coordinates in a neighborhood of O and in a neighborhood of O, denoted {q,p) = (<7i, Q2, Pi, P2) 
and (q,p) = (qi, fa,Pi,P2), having the origin at O and O respectively. It is clear that the change (13j) relating them is 
given by 

(Qu fa) = X(H, Q2) = (qi - 27r,g 2 ), (Pi>P2) = (pi,Pa)< (58) 

since D\(q)~ T = Id. 

Remark. Due to the periodicity of the Hamiltonian, the corresponding Hamiltonian equations write identically in 
either coordinates (q,p) or (q,p). The invariant manifolds W u ' s and W u ' s are related by the periodicity and, therefore, 
their generating functions S u,s and S u,s are the same functions (note however, that the functions S u ' s depend on q 
whilst S u,s depend on q). 



We assume the hypotheses (HI— H4) stated in Section [O] with some minor modifications in order to adapt them 
to the case Q = T 2 being studied here: 

(H2') the orbit 7 satisfies (72 = 0, with qi increasing along the orbit from to 27r; 

(H4') the generating function S u (qi,(7 2 ) of the invariant manifold W u of O is defined in a M 2 -neighborhood of the 
segment < qi < n of type U = (—a, n + a) x (—5, S), for some a > and S > 0; and conversely, the generating 
function S s (qi,fa) of the invariant manifold W s of O is defined in a M 2 -neighborhood of the segment — n < q\ < 
of type U = (— 7r — a, a) x (— 6, S). 



Applying the change (58), we can write the manifold W s in the coordinates (q,p). As in (|6|, it has a generating 
function, given by S s = S s o x + const = S s o \ + const, that is, 



5 s ((? 1 ,(72) = 5 s ((7 1 -27r,g 2 )+(7. 



(59) 



In other words, the generating function S s (q) is the function S s (q) expressed in the coordinates q and shifted by an 
additive constant a, which will be chosen below. 

By hypothesis (H4'), we see that S u (q) and S s (q) can be compared in their common domain U n X~H^) = 
(it — a. 7r + a) x (—5, 5), which is a R 2 -neighborhood of q\ = it. As we see below, our choice of this point to check the 
transversality allows us to take advantage of the reversibility properties of the Hamiltonian. 



The generating functions along the loop. As in the general setting described in Section |2.2[ the loop 7 is 
described by the coefficients S and Si, of orders and 1, in the expansion in q 2 of the generating function of any of 



the invariant manifolds containing it. More precisely, the loop 7 is given by the equations (24-25), which correspond to 



the coefficients in the expansions ( 19 ) and ( 23 ). Such coefficients are, initially, functions of q\ defined in a finite interval. 



Nevertheless, using that they must coincide for the unstable and stable manifolds, together with the 27r-periodicity 
in q\ of the Hamiltonian, we are going to show that such coefficients are 47r-periodic functions of q\. 



Combining (24 25) with the identity (59), we have 



7 : 



Pi 

P2 



S' (qi) 
Si(gi) 



(<So)'(<?i) - v^o 
Si u (gi) = St(qi) 



(S s )'(qi) 



(S s )'(qi 



2tt), 



(60) 
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As a parameterization of 7, this is valid for < q% < 2tt, but the functions in the equalities can be extended up to a 
neighborhood of < qi < 2tt. 



In fact, it is easy to see from the above equations (60) and the reversibility relations of Proposition 10 with 
ri = r% = 1 — this is the trivial 7?.-reversibility, see the remark after Lemma [9] — , that S' and Si can be extended as 
27r-antiperiodic functions: 

S'o(qi) = (Sg)'(«i - 2tt) = -(S u )'(<?i - 2n) = -S ( qi - 2tt), 
Si( qi ) = Sl( qi - 2tt) = -Sf( qi - 2tt) = - 2tt). 

Thus, the functions S' (qi) and Si (51) can be defined for all <7i G R, and are 27r-antiperiodic, and hence they are 
47r-periodic and have zero average. It turns out that the coefficient So(<Zi) (without derivative) is also a 4-7r-periodic 
function. Indeed, one can easily check that a primitive of any 27r-antiperiodic function is 47r-periodic. 

We point out that this period 4tt is due to the fact that the functions S' a (qi) and Si(gi) parameterize two loops 
of O. Indeed, the loop 7 is obtained for < qi < 2ir and satisfies pi > 0, and another loop is obtained for 2ir < qi < Air 
and satisfies p\ < (see an example in Figures 2(a) and 2(b) ). 



Now it is clear, by integrating from to 2tt the first equation in (60), that we can write 

Sofa) ■= S u (-7i) - S s ( qi ) = S s ( qi - 2vr) + a, 



where a is the constant in ( 59 ) , given by 

/.2tt 



o=\ pidgi = S (2^) - S o (0) = S u (2tt) = -S^(-2tt), 
Jo 



(61) 



and we have taken into account our choice Sq' s (0) = in (19). 



The generating functions of the invariant manifolds. Let us return to the invariant manifolds W u and W s , 
described in the coordinates (q,p) by their generating functions S u (91,^2) and S s (qi,q2). We know from the previous 
paragraph the relations concerning the coefficients of orders and 1 in their expansions in qi , and now we are interested 
in the coefficients of order 2, which are used in order to study the transversality. 



Expanding in qi the identity (59), the coefficients of order 2 satisfy the equality 

f s ( (Zl )=T s ( (Zl -27r). 




Figure 2: (a) Graph of the 2n -periodic potential Vo (qi) = cos qi — sin qi + cos 2qi — | sin 2qi — 2. (b) Projection onto the plane 
(qi,pi) of the loops, given by the graph pi — Sb(qi) and assuming P{qi) = 1 (see Proposition^^b)): the loop 7 has qi increasing 
from to 2-7T, and a second loop obtained from the periodicity has q\ decreasing from <in to 2tt. 
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By hypothesis (H4'), the functions T u and T s are defined for q x G (—a, tt + a) and qi E (tt — a, 2tt + a) respectively. 
Then, they can be compared for q\ G (tt — a. tt + a), which makes a difference with the fact that the coefficients of 



orders and 1 (giving only the loop in (60)) are functions defined for all q\ G 



Reversibility relations and transversality. Finally, we assume the reversibility of the Hamiltonian in the sense 
sta ted in Section 2.5 i.e. with respect to an involution 1Z of the form (38 1 with the matrix R = diag[r!,r 2 ] given 
in (39). Since we want to compare W u with W s along 7, we shall fix r\ = — f, for this will be the case of interest: 
the loop is mapped into itself by the reversibility, TZ^f = 7, since it reverses the sign of q±, but keeps the sign of p±. In 
other words, for n = — 1 the restriction of 1Z to the plane (q\,pi) is a symmetry with respect to q\ = (or cquivalcntly 
qi = tt). 



Proposition 13 Under the hypothesis of this section, if the Hamiltonian |Ip is IZ-reversible with 7*1 = 
r 2 = H), then the following equality holds: 

S s (2ir- q\,r 2 q 2 ) = -S u (qi,q 2 ) + a, (qi,q 2 ) € (tt - a, tt + a) x (-6,5) , 



— f (and 
(62) 



where a is the constant in (59). Besides, its coefficients satisfy 



S (2TT- qi ) = -Sofa) + a, qi€R, 
Si(2rr-gi) = -r 2 Si(qi), qi G K, 

f s (2TT~ qi ) = -T"( qi ), qi e(TT-a,TT + a). 



(63) 



Proof. We use the reversibility relations given by Proposition fO with r\ = — 1. The first equality (62 1 follows 



from (59) and the first relation in Proposition fO 



S s (2tt - qi,r 2 q 2 ) = S s (-q 1 ,r 2 q 2 ) + a = -S u (qi,q 2 ) + a. 



Expanding this equality in q 2 , we get the remaining three equalities (63) 



However, to ensure the validity for q± ^ (tt — a,TT + a) of the first two equalities <|63]) , we should deal with the 
whole loop. For a given point (qi,0,pi,p 2 ) G 7, we have in (60) the equations p\ — S$(qi) and p 2 = Si(qi). Since 1Z 
takes this point into another point of 7, we also have p\ = S' q (2tt — qi) and —r 2 p 2 = Si(2tt — qi). Then, the functions 
So and Si satisfy the identities S' q (2tt — qi) = S' (qi) and Si(27r — qi) — —r 2 Si(qi). The second one corresponds to 
our statement, and integrating the first one we obtain So^tt — 9i) = — So(<Zi) + const, where for q ± = we see that 
the constant is a, in agreement with (6f ). □ 



As a direct consequence of this proposition, we see that in this case Q = T 2 the transversality condition (Tsl) 



becomes very simple. Using the third equality (63) at q\ = tt, we get that the necessary and sufficient condition for 
the transversality of the invariant manifolds W u and W s along the loop 7, can be expressed as follows: 



T u (tt) y£ 0. 



3.2 Example: two identical connected pendula 

As an example with the configuration manifold Q = T 2 , we consider the system of two identical pendula connected 
by an interacting potential. This model was considered in |GS95j . where the transversality of the invariant manifolds 
along a loop was established by using the variational equations around it (for other properties of this model, see 
for instance |Arn891 §23D]). Here, our aim is to apply the method developed in this paper to recover the result of 
transversality in a somewhat more general model. 

Initially, in symplectic coordinates (£, 77) = £ 2 > Vis V2) <= T 2 x W 2 we consider the Hamiltonian 

Hfori) = Uv? + Vi) + (cosCi - 1) + (cos 6 - 1) + /(&)(! - cos(6 - &)). (64) 
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where the interacting term is given by a 27r-periodic even function /(£i), such that < /(0) < 1/2. Under this 
inequality, the origin is a hyperbolic equilibrium point with (different) Lyapunov exponents ±1, ±6, where we define 
b = y/1 - 2/(0), < b < 1. The case of a constant function /, considered in [GS95] , corresponds to a linear spring 
connecting the two pendula. We point out that the term 1 — cos (£2 — £1) allows us to keep the periodicity in £1, but 
could be replaced for instance by (£2 — £i) 2 /2. 



It is clear that, if we assume / = in (64), the system is separable and integrable: it consists of the product of 



two identical pendula with coordinates {£j,T]j), j = 1, 2. We have, for each pendulum, the (positive) separatrix 



£j = 4 arctan e* , rjj = 2 sin ■ 3 



(65) 



2 cosh t 

Combining the two separatrices, with a free initial condition for one of them, say the first one, we have the following 



1-parameter family of loops indexed by s £ M. (see Figure 3(a) ) 

7s : 



£i(i, s) = 4 arctan e s+t , £a(£) = 4 arctan e*, 

2 2 
Vi(t,s) = — — —r , fj 2 (t) 



(66) 



cosh(s + 1) 



cosh t 



If we consider / ^= in (64), in general the family of loops is destroyed, but the concrete loop 7 = 70 remains 



unchanged, as well as the dynamics along it. Our aim is to provide a sufficient condition on the function /, ensuring 
that the invariant manifolds intersect transversely along this loop. 



In order to apply the results of Section [2] we perform the symplectic change of coordinates 

A>P = I 'V I ■ I " I =B l p(q)- T p : 



6 



qi + <?2 



m 
m 



Pi ~P2 
P2 



(67) 



where it is clear that £ = </?(<?) is a well-defined change on T 2 . In the new coordinates (g,p), the expression of the 
Hamiltonian H takes the form ([I]), with 



B(q) = 



1 -1 
-1 2 



It is easy to see from Lemma [9] that H is 7£-reversible, with n 
the loop 7 lies on the line q 2 = 0: 



, V(q) = (cosgi - 1) + (cos(<?i + q 2 ) - 1) + /(?i)(l - cosg 2 )- 

1 (recall that / is an even function). Now, 



«?(*) = 4 arctan e*, ^(t) = 0, p\(t) = p° 2 (t) = 



cosh t 



We know from ( 24 ) that the loop is also given by 



<72=0, pi = S' (qi) =4shiy, p 2 = Si(qi) = 2sin y . 



Notice that these are 47r-periodic functions, in agreement with the considerations of Section 3.1 



Remark. As an alternative to (67), we might have chosen as £ = ip(q) a rotation taking the line £1 = £2 into 
q 2 = 0. This change would lead to a diagonal matrix B(q) and, although it is not a one-to-one change on T 2 , some 
arrangements can be done. However, we chose the change (67) since it is a particular case of the change considered 



for the example of Section 4.2 



Now, we can write down the Riccati equation ( |34[ ) allowing us to study the transversality. Developing V(q) in q 2 , 
we get: 

V (qi) = 2(cosgi - 1), Vi(gx) = - singi, Y(qt) = cosgi - f(q t ). 
Then, the Riccati equation for the unstable manifold W u is 



2 sin I (T u )' _ 2 cos I T u + 2(T U ) 2 = - (sin 2 | + /(«&)) , T u (0) = - 



(68) 
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where the initial condition comes from Lemma [7J We give below a sufficient condition on /, ensuring that the solution 
of (68) can be extended up to qi = it, and satisfies the transversality condition T u (ir) ^ 0. 



We provide another useful form for the Riccati equation (68 1 by carrying out the change of variable x = — cos(q±/2), 



with a further modification in order to remove the linear term. Writing T u (q\) = (T n (x) — x)/2 and f(qi) — f(x), we 
get the equation 

(l-* 2 )— + (T u ) 2 = -2(l-* 2 ) + (l-2/>)), T u (-1) = 6, (69) 
and the transversality condition becomes T u (0) ^ 0. 



As in Section 



2.6 



second-order linear differential equation 



- 1 — x 2 dy 

the Riccati equation (68) can be transformed, through the change T u (x) = — — — — , into the 



y(x) dx' 



d^y 
dx 2 



2x 



dy 
dx 



1 - 2 / (x) , 



(70) 



We are going to establish the transversality of the invariant manifolds along the loop 7 by checking the transversality 



condition, for a wide class of functions / in (64). We start with the simplest case of a constant function, already 
considered in GS95], 



Lemma 14 If we consider 1 — 2f(x) = b constant, with < b < 1, then the solution of (69) is 



n{ X )=b- 



l~x l 



(71) 



which can be extended up to the interval x G [— 1, b), and satisfies T^(0) = b — 1/b < 0. 



Proof. In this case, the second-order equation ( 70 ) is a generalized Legendre equation, for which a fundamental system 
of solutions is given by 



2/1(2:) = (b + x) 



1 - x 



b/2 



y 2 (x) = (b-x) 



1 + x 
1-x 



6/2 



= Vi{-x). 



We get, from the second one, the solution of (69) as given in (71). This solution obviously satisfies the transversality 



condition, provided < b < 1 (of course, the separable case b = 1 has to be excluded). 



□ 



Proposition 15 If the function f satisfies < f(x) < 1/2 for all x G [—1, 0], then the solution T u (x) of \ 69) can be 
extended up to the interval x G [—1,0], and satisfies T u (0) < 0. 



Proof. The argument is similar to the one given for the example of Section 2.6 With a given f(x), we write (69) in 
df u 

the normalized form — — = F(x, T u ), where the function F is defined for (x, T u ) e (—1, 0] x R. Let c and d be positive 



d.j- 



constants such that < c 2 < 1 — 2f(x) < d 2 < 1 for any x. Replacing 1 — 2f(x) by c and d in (69 ), we define functions 
F c (x,T u ) and Fd(x, T u ), corresponding to respective differential equations of the type considered in Lemma 



14 



We 

know that their associated solutions T" and T£ are both defined for x G [— 1, c) since c < d. Combining the initial 
conditions T c u (-1) = c, T u (-1) = b, T^(-l) = d, that satisfy c < b < d, with the fact that the inequality F c < F < F d 
is true for x > —1, we see that the solution T u (x) is confined between T^{x) and T^(x), and hence it can also be 
extended up to the whole interval [—1, c), and satisfies T u (0) < T^(0) < 0. □ 



Remark. If the function / in ( 64 ) is a trigonometric polynomial, then the change of variable between q\ and x leads 



to a second-order linear differential equation (70), whose coefficients are rational functions of x; hence, as in Section 2.6 



one can apply Kovacic's algorithm to see whether (70) is integrable (in the Liouville sense) or not and, when it is, 



obtain explicit solutions. In such a case, the transversality condition could be checked directly. Nevertheless, applying 
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Proposition [15] we can check the transversality in many cases, even when it is not possible to solve ( 70 ) explicitly. 
As a concrete example, choosing = (2 — cos£i)/8 in (64), we have fix) = (3 — 2x 2 )/8, and applying Kovacic's 



algorithm (as implemented in Maple) for this case it turns out that the equation (70 1 cannot be solved explicitly but, 



on the other hand, without solving it we see from Proposition 15 that, for this example, the transversality condition 
is fulfilled. 



4 The case of a perturbed Hamiltonian 

4.1 Mel'nikov integrals 

Now, we study the invariant manifolds in a Hamiltonian depending on a small perturbation parameter e, 

H s (q,p) = H(q,p)+eH*(q,p). (72) 
As in Section [3] we assume that the configuration manifold is Q = T 2 . For the unperturbed case e — 0, we assume 



that the Hamiltonian Ho = H is in the situation described in Section 3.1 In particular, it has a hyberbolic point O 
at the origin, with a loop 7 contained in q 2 = 0. 

Concerning the invariant manifolds in a neighborhood of the loop, we are going to consider two different situations: 

(A) the unperturbed invariant manifolds W u - S intersect transversely along the loop 7; 

(B) the unperturbed invariant manifolds W u ' s coincide along the loop 7. 

We stress that loops of both types (A) and (B) may coexist in the same Hamiltonian H (and each of them can be 



moved to q 2 = by a specific symplectic change), as illustrated by the example of Section 4.2 



Notice that a loop 7 of type (A) is isolated in the sense that, if an orbit close to 7 is a loop, it is not contained in 
a small neighborhood of 7. 

On the other hand, if 7 is a loop of type (B) it belongs to a 1- parameter family of loops, filling the 2-dimensional 
homoclinic manifold or separatrix W := W u = W s . 

For e ^ 0, in both cases (A) and (B) the loop 7 will not survive in general. Our aim is to study whether, for e small 
enough, there exists a perturbed loop r y £ at distance 0(e) from 7, and the perturbed manifolds W" and W| intersect 
transversely along -f e , inside the energy level containing them. We show that, in the case (A), this problem is quite 
simple and can be solved by applying the implicit function theorem on a transverse section q\ — const. Instead, in the 
case (B) we have a problem of splitting of separatrices, which must be re-scaled in order to apply the implicit function 
theorem, leading to a necessary and sufficient condition which can be expressed in terms of Mel'nikov integrals. In 
other words, we have for (A) and (B) a regular and a singular perturbation problem respectively. 

It is clear that for e small enough we have a perturbed hyperbolic equilibrium point O e inherited from O, whose 
invariant manifolds W"' s can be expressed in terms of perturbed generating functions: p = V5"(g) and p — VS £ (q) 



(see Section 2.2 for the notations). We point out that the identity (59 1 also applies here to the perturbed case: 



Sl(q 1 ,q 2 ) = S s e (q 1 -27r,q 2 ) + a £ . 

Let us consider the expansion in s of the generating functions, 

S»( q )=S"( q )+eS^q) + 0(e 2 ), S s e (q) = S s (q) + eSt{q) + 0{e 2 ). (73) 



On the other hand, we know from the ideas introduced in Section |2.2| that the transversality of the perturbed manifolds, 
inside the energy level containing them, can be studied from the expansion in q 2 , up to order 2, of the generating 



functions (73). For the unstable manifold, we write 



S»(q) = Sl e ( qi ) + SIM q 2 + \T?{ qi ) q 2 + 0(q 2 3 ) 
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with 

Sf^) = S.iq^+eSlM+Oie 2 ), j = 0,1, 
T?(q 1 ) = T»(q 1 )+eT:( qi )+G(e 2 ), 

and analogously for S £ (q). Notice that, for e — 0, we have written Sj(qi) :— Sj(qi) = Sj(qi), j = 0, 1, since these 
terms coincide due to the existence of the unperturbed loop 7 (recall the definitions (25)). 

Following a standard terminology (see for instance IX illl)] ). we define the splitting potential as the difference 
between the generating functions for the unstable and stable manifolds: 

C e (q)=AS s (q):=S2(q)-S B e (q) (74) 

(for any functions / u and / s , we are using the notation A/ = f u — f s ). We point out that the name 'splitting potential' 
is normally used in the context of splitting of separatrices, i.e. our case (B), but we are considering here a somewhat 
more general situation. Expanding C £ in q 2l we write 

C e {q) = AS ,s(<7i) + &SxM 92 + §AT £ ( ?1 ) di + 0(qi). 

If we consider the expansion C e = C + e£* + 0(e 2 ), we have 

C( q ) = ^AT( qi )qi+0( q i), 

C*(q) = AS ,*(<?i) + A5i,,(9i)ga + h AT M li + 0(qi), 
where we used that AS'o(<Zi) = A5i(gi) = 0, by the existence of the unperturbed loop 7. 

Now, we can rewrite the two hypotheses on the loop, in terms of the unperturbed generating functions. Recalling 
the transversality condition Q for the case (A), we have: 

(A) AT( qi ) f for any ?1 ; 

(B) AS(q) = for any q=(qi,q 2 ). 



Theorem 16 

(a) If j is a loop of type (A), then for e small enough there exists a perturbed loop j e biasymptotic to O e , and the 
invariant manifolds W"' s intersect transversely along 7 e . 

(b) If j is a loop of type (B) and the conditions 

°P- (tt, 0) = ASr,* (tt) = 0, ^(tt, 0) = AT, (tt) ^ (75) 
d q2 dq 2 z 

are fulfilled, then for e 7^ small enough there exists a perturbed loop 7 e biasymptotic to O e , and the invariant 
manifolds W" ,s intersect transversely along "f e . 



Proof. As in Section |2.2[ the energy level Af s containing both perturbed invariant manifolds is given by the equation 
H £ — H £ {O e ) = 0. Since for e = we have = 5i 7^ on 7, we see from the implicit function that, near 7 and for 
e small enough, we can parameterize Af e as p\ = g e (qi, q 2 ,p 2 )- I n the coordinates (c[i, q2 ,p 2 ) of Af s , the intersections 
between the invariant manifolds W" and W| are given by the equation 




It is enough to consider the intersections in a transverse section, say q\ = tt, and solve the equation for q 2 . Expanding 
in q 2 , we see that the intersections are given by the solutions of 

j-^(n, q 2 ) = A5i, e (7r) + AT e (7r) q 2 + 0(q 2 2 ) = 0. (76) 
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Now, expanding in e and recalling that ASi — 0, the equation becomes 

(AT(tt) q 2 + 0(qi)) + e (A5 X ,, (tt) + AT* (tt) g 2 + 0(g 2 2 )) + 0(e 2 ) = 0. (77) 

Therefore, in the case (A) there is a solution q 2 = 0(e), and hence the invariant manifolds intersect transversely along 
a perturbed loop j £ , for e small enough. 

In the case (B), the fact that C(q) = AS(q) = for any q implies that — — (^,q 2 ) = AT(tt) q 2 + 0(q 2 ) — 0, and 



hence the expansion (77 1 begins with terms of order 1 in e. Then, dividing by e the whole equation we obtain: 

AS h *(n) + AT,(tt) q 2 + 0( 9 2 ) + 0(e) = 0. 



Since we assumed in (75| that A5i.*(7r) = and AT*(7r) 7^ 0, we see that there is a solution q 2 = 0(e), and hence the 



invariant manifolds intersect transversely along a perturbed loop 7 e , for £ ^ small enough. □ 

Remark. In the case (B), it is easy to see from the existence of the perturbed loop 7 E that A5i i *(gi) = for any q\. 

Indeed, parameterizing 7 e by qi, we write q 2 = f £ (qi). It is satisfied that E (qi, feili)) — 0. Expanding this equation 

oq 2 

in e and using that for e = we have C — and f(qi) — 0, we obtain: A5i.*(gi) = (gi,0) = 0. 

oq 2 

Remark. We point out that, for a loop of type (B), it is not enough to impose a condition like ([8]), and we need 
as an additional condition that the coefficients ^,(tt) and S\ ^{ir) coincide. Otherwise, the existence of a perturbed 
loop 7 £ close to 7 could not be established in general. 

The remaining part of this section concerns the case (B), and our aim is to provide an explicit condition allowing 
us to check ( |75| in a concrete example. For e = 0, the separatrix W is a graph p — VS(q), with the generating 
function S(q) := S u (q) — S B (q), and the inner dynamics on VV is given by the ordinary differential equation 

q = B(q)VS(q). (78) 



For a given q = (qi,q 2 ), denoting q(t,q) the solution of (78 1 that satisfies the initial condition q(0,q) — q, we have a 
loop 

x{t,q) = {q(t,q),P(t,q)), with p{t, q) = VS(q(t, q)), (79) 
Of course, it will be enough to consider the initial condition in a direction transverse to q 2 = 0: 

q = k(s) = (ki(s),k 2 (s)), se(-a,a), with k(0) = (7t,0) and k' 2 (0) ^ 0. (80) 

In this way, we have a 1-parameter family of loops indexed by s, that we denote j B , given by q(t, s) := q(t, k(s)). The 
particular loop contained in q 2 = is 70 = 7, and we write q (t) = (q^(t),0) — q(t,0). Notice that (qi,q 2 ) — q(t,s), 
(t, s)Glx (—a, a), is a change of parameters for the separatrix VV, in a neighborhood of 7. 



As a possible approach, we could see the functions T"(gi) and T^(q±) as solutions of Riccati equations as in ( 34 ) , 
with initial conditions at qi = and q\ = 2tt respectively. Developing such equations in e and taking the terms 
of order 1, the functions T"(gi) and T%(qi) become solutions of linear differential equations, and hence they can be 
expressed in terms of integrals. Comparing such integrals, one would obtain a condition for the transversality which 
could be checked. 

Nevertheless, similar integrals of Mel'nikov type can be obtained, in a more direct way, by developing in e the 



Hamilton- Jacobi equation satisfied by the perturbed generating functions ( 73 ) and comparing the unstable and stable 
ones. An analogous approach was previously used in [LMS03 (see also |Sau01| ). for the case of invariant manifolds 
of tori associated to a simple resonance, i.e. with a 1-d.o.f. hyperbolic part. We are going to show that a first order 



approximation for e small of the splitting potential C e defined in (74 1 is given by the Mel'nikov potential, defined as 
the integral 

/•oo 

L(q) = - [H, {x(t,q))- H, (O)] dt, (81) 
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absolutely convergent due to the fact that x(t,q) tends to O, with exponential bounds, as i ^ ±00. Essentially, 
the Mel'nikov potential is defined as the integral of the perturbation on the unperturbed loops. See also |Tre94] 
for another definition of a Mel'nikov integral along the loops contained in a separatrix, and |DR97j for the case of a 
hyperbolic fixed point of a symplectic map. 

Theorem 17 The splitting potential is given at first order in e by the Mel'nikov potential: for some constant c, 

£ E (q)=e(L(q) + c)+0(e 2 ). 



Proof. Expanding in e the generating functions in (74), it is enough to show that the difference of their first order 
terms is given by the Mel'nikov potential: 



S»{q) - S S M = L(q) + c. 



(82) 



As in the proof of Theorem |6j we use that the perturbed generating functions satisfy the Hamilton- Jacobi equation, 
but now we expand it in e. We start with the generating function associated to the unstable manifold: for any q, we 
have 

H e (q,WS i ;(q))=n(q)+en:(q) + 0(e 2 )^eH,(0)+0(e 2 ) 

(to obtain the constant at the right hand side, we expand H £ (0 £ ) in e and we use that H(0) = and VH(0) = 0). 
The terms of orders and 1 are 

H(q)=H(q,VS(q)), 

Kiq) = (^(q,VS(q)),VS»(q)\+H*(q,VS(q)). 

Since for any q we have H*{q) = H*(0), we deduce from that VS 1 " (q) satisfies the following linear partial differential 
equation: 

(B(q)VS(q),VS«(q)) = -{H*(q,VS(q))-H40)}. (83) 
In the solution of this equation, the characteristic curves play an essential role. Since they are given by the ordinary 



differential equation ( 78 ) , the characteristic curves are projections of the unperturbed loops ( 79 1 onto the configuration 
space: q = q{t, q). 

What follows is a standard argument for the Poincare-Mel'nikov method. Using that Sl l (q(t, q)) tends to 5"(0, 0) 
as t — > — 00 for any q, we see that 



,S". l Uy) >"'((>.(>) / ±[S?(q(t,q))]dt 



-OO 





(VS3f($(t,«)),$(i,g)) dt 



[H*(x(t,q))-H*(0)} dt, 



where we have used (78) together with (83). Proceeding analogously with the generating function associated to the 

/>oo 

Sl(2w,0) - Sl(q) = - / [H.(£(t,q))-H m (0)]dt, 



stable manifold, we have 



and adding the two expressions we get the identity (82), with c = 5" (0,0) — 5^(2^,0) 



□ 



We see from this result that, in terms of the Mel'nikov potential, the transversality condition (75) becomes 

d 2 L 



dqi 



(7r,0)=A5i-(7r) = 0, 



dq. 



2 (7r,0) = AT t (7r)^0. 



(84) 



Now, we show in the following easy lemma that the Mel'nikov potential L is constant along each loop or, in other 



words, it is a first integral of the inner dynamics on the separatrix, given by (78) 
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Lemma 18 The value of L(q(r,q)) does not depend on r. 



Proof. It is enough to carry out the change of variable t H> t — t in the integral (81 1, together with using the identity 
$(t,q(T,q)) = x(t + T,q). □ 



As a consequence, the splitting potential can be studied in terms of the parameter s introduced in (80 1. With this 
in mind, and denoting x(t, s) := x(t, k(s)), we define the reduced Mel'nikov potential as 



L(s) := L(k(s)) 



[H m {x(t,s))-H,(0)]dt, sg(-o,o) 



(85) 



We show in the next result that the transversality condition can also be written in terms of L(s) (for an illustration, 
see the example of Section 4.2). 



Proposition 19 If the reduced Mel'nikov potential L(s) has a nondegenerate critical point at s = 0, then for £ ^ 
small enough there exists a perturbed loop j e biasymptotic to O e , and the invariant manifolds W"' s intersect transversely 
along r ) e . 



Proof. We first show that 

Z'(o) = o 

and, when this is satisfied, we also show that 

Z"(0) ^ 



dq 2 



M) = o 



d 2 L 

dq£ 



(tt,0)^0. 



Then, the transversality condition (84) is fulfilled when the reduced potential L(s) has a nondegenerate critical point 
at s = 0. 



To show the two equivalences, we use the identity 

L(q(t,s))=L(s), 



(86) 



which comes from Lemma 18 and the fact that q(0, s) = k(s). The first derivative with respect to s of the identity (86) 
provides 

-ms))^ { t,s) + - m , S ))^(t,s)^L'(s), 

Considering s — and using that L is constant along q 2 — 0, we get the equality ^—(q(t, 0)) 0) = L'(0). Since 

oqi as 

by (80) we have -^(0, 0) = Ko(0) 7^ 0, we obtain the first equivalence. Now, assuming that L'(0) =0we deduce that 
os 

d 2 L , d 2 L 



(q(t,0)) = at least in a neighborhood of t = 0. Thus, we have — — 2(1", 0) 

oq-y 



8L 
dq 2 

second derivative with respect to s of (86), at s = t = 0, provides 

8 2 l 



dqi dq : 



(7T, 0) = 0, and hence the 



dq. 



which implies the second equivalence. 



2 (tt,0) 4(0) 2 = L"(0), 



dC, 



□ 



Remark. If one defines the splitting function as M. e (q) = — — (q), this provides a measure of the distance between 

oq2 

the invariant manifolds in the p2-dircction. At first order, this function can be approximated in terms of the Mel'nikov 

dL ~ ~ 

function M(q) = — — (q), or the reduced one M{s) = L'(s), which has to have a simple zero at s = 0, as the condition 

for the transversality of the invariant manifolds. 
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To end this section, we revisit the above results assuming for the perturbed Hamiltonian a suitable type of re- 
versibility (see Section 2.5 1. 



Proposition 20 If the Hamiltonian H e is 7?.-reversible with r\ = r 2 = — 1, then £ £ (ir, q 2 ) and L(tt, q 2 ) are even 
functions in q 2 . 



Proof. We see from the considerations of Section [37lj in particular Proposition 13 that the generating functions of 



the invariant manifolds ,s are related by S s e {2ir — qi, — q 2 ) = — S£(qi, (72) + o e . Then, we obtain 

£ £ 0r,<z 2 ) = S"(ir,q 2 ) + Sz(ir,-q 2 ) + cr e , 
an even function. This applies also to the Mel'nikov potential L(ir,q 2 ), though one could also check that this is an 



even function directly from its definition (81 1. □ 



Clearly, the fact that C e (ir, 92) is even in q 2 implies that the equation ( 76 ) always has q 2 = as a solution, implying 
the existence of a perturbed loop "f e (whose (72-coordinate coincides with 7 at qi — n but, in general, we will have 
7 £ 7^ 7). Nevertheless, an additional condition has to be imposed in order to ensure the transversality of the invariant 
manifolds along the loop j e . 



Additionally, if we choose q — n(s) in (80) in such a way that the family of unperturbed loops satisfies 7_ s = 7Vy s , 
then the reduced Mel'nikov potential L(s) is an even function in s. Then, it always has a critical point at s = 0, and 
we only have to check its nondegeneracy, L"(0) ^ 0, in order to ensure transversality. 



4.2 Example: two different weakly connected pendula 

In symplectic coordinates (£,17) G T 2 x R 2 we consider the Hamiltonian 

He(t n) = + tfc 2 ) + (cos£i - 1) + A 2 (cos£ 2 - 1) + e(l - cos(& - £1)), 
where £ is a small parameter, and A > 1 is a fixed value (recall that the case A = 1 has already been considered in 



Section 3.2 1. 



For e = the system is separable, and consists of two pendula with Lyapunov exponents ±1 and ±A, generalizing 
the separable case of Section |3.2| On the region 771, 7/2 > (the other three ones are symmetric), it has a 1-parameter 
family of loops 7,, s 6 E, plus two special loops 71 and 72 with a different topological behavior. The two special loops 
are given by the separatrix of first/second pendulum, with the equilibrium point of the second/first pendulum: 

71 : £1 (t) = 4 arctan e* , r?i (t) = , & (i) = m (t) = 0, 

At 2A 



72 : £1 (t) = m (t) = 0, 6 (t) = 4 arctan e At , r? 2 (t) 



cosh Xt 



and the 1-parameter family is given by the separatrices of the two pendula, with a free initial condition in one of them 
(see Figure [3(b)] ) : 

% : £1 = s) = 4 arctan e t_s , £ 2 = £2 (t) = 4 arctan e At , 

2 2A 

771 = fjx(t, s) = T7- , 772 = f} 2 (t) 



cosh(7 — s) ' cosh Xt ' 

It is not hard to see that the loops 71 and 72 are of type (A). Indeed, the unperturbed unstable/stable invariant 
manifolds of the loop 71 are given by the separatrix of the first pendulum, times the local unstable/stable curves of 
the equilibrium point of the second pendulum: recalling ( |65| , their equations are 771 = 2sin(£i/2), 772 = ±2Asin(£2/2) 
(with the signs +/— for the unstable/stable manifold respectively), and it is clear that they intersect transversely. 



Then, we see as a consequence of Theorem 16 a) that the equilibrium point O e = O has a perturbed loop 71 E whose 



invariant manifolds intersect transversely along it for e small enough (from the Hamiltonian equations, one sees that 



31 



for £ ^ there is no orbit contained in £ 2 = 0, and hence 71. E 7^ 71). The same considerations are valid for the loop 
72, since we did not use the fact that A > 1, and both loops are completely analogous. 

Now, we study the loop 7 = 70, corrsponding to s = in the 1-parameter family. The existence of this family 
implies that this is a loop of type (B), whose separatrix W is the manifold defined by the equations r\\ = 2sin(£i/2), 
rj2 = 2Asin(^ 2 /2), < £1,62 < ^ 7r - The manifold W contains the whole family of loops j s , s € R. 

In order to introduce new coordinates (q,p) such that the loop 7 is contained in q 2 = 0, we take into account that 
the ^-projection of the loop 7 is the graph 

£2 = h(£i) '■= 4arctan ^tan A ^- 

This function h : T — > T is of class C r , where r is the integer part of A, and h is analytic if A is integer. Notice that 
we can consider h as an odd function, since the equality h{2ir — £1) = 2-k — is fulfilled. 

Thus, we consider the change of coordinates 

\t) =iP{q)= \h^) + q ^ 



(notice that x = ip(q) is a well-defined change on T 2 ). In the new coordinates (q,p) the Hamiltonian takes the form ( 72 ) 
where, for the unperturbed part H , we have 



1 -h'{q x ) 
^ -h'( qi ) l+h'( qi ) 

V(q) = (cosgi - 1) + X 2 (cos(h( qi ) + q 2 ) - 1), 
and the perturbation eH* is given by 

fl,(g) = l-cos(fc(? 1 )-?i + gjj). (89) 
Let us write down, for £ = 0, the q-coordinates of the unperturbed loops, 

j s : g x = |" x = 4arctane*~ s , q 2 = £2 - Mfi) = 4arctane At - 4arctane A(t_s \ (90) 



and it is clear that the loop 7 = 70 is now contained in q 2 = 0. As in (80), the initial conditions in (90) are given by 
k(s) = (4arctane _s , ir — 4 arctane _As ), a direction transverse to q 2 = for s = 0. 
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As we easily check, the perturbed Hamiltonian H E in (87 89 1 is 7?.-reversible with r\ = r 2 = —1 (using that the 



functions h and h! are odd and even respectively). Then, we see from Proposition 20 that there exists a perturbed 
loop 7 e close to 7. 

To study the transversality of the perturbed invariant manifolds along we use the reduced Mel'nikov potential: 

/oo />oo 
HM dt=- [1 - cos(/i(g 1 ) - q x + q 2 )} dt 
-OO J —OO 

/OO 
[1 - cos(6 - &)] dt, (91) 
-co 



where the expressions for £1 = £i(t, s) and £2 = &(t) have been introduced in (66 1 



The following result implies, according to Proposition |19[ that the invariant manifolds intersect transversely along 
the perturbed loop 7 e for e 7^ small enough, at least for the interval of A considered. 



Proposition 21 Assuming that 1 < A < Ao — 3.68078, the reduced Mel'nikov potential (91) has a nondegenerate 
critical point at s = 0. 



Proof. Differentiating (91) with respect to s, we obtain: 



L'(M=/ sin(6-^)^dt, 



where we have 



ds 



L"{s) = 
-2sin(£i/2) and 



^1 
ds 2 



cos(6 - £1) 



sin £1 . 



ds 



sin(£ 2 - £1) 



ds 2 



(92) 



dt, 



Due to the 7?.-reversibility and the fact that the family of unperturbed loops satisfies 7-s = T^-ls, it turns out 
that L(s) is an even function, and hence L'(0) = (this is a consequence of Proposition 20 or may also be checked 
directly) . 

Now, let us see that £"(0) < 0, at least for a large interval of values of the parameter A. We are going to show 
that for s = the function inside the integral in ( 92 1 , 



f(t) = 4cos(6 - sin 2 I - sin(£ 2 - £1) sin fx, 



is positive for any t € R. Notice that for s = we have £1 = 4arctane* and £2 = = 4arctane 



A/ 




Figure 4: 27ie reduced Mel'nikov potential L(s) for A = 1. 
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For any t > 0, the difference £2 — £1 is positive, and reaches its maximum value S> when t is such that cosh At = 
A cosh t. Using that £2 = 1S increasing in A for t > 0, then we see that S> is also increasing from to ir as A goes 

from 1 to 00. Therefore, there exists a value Ao such that S> = ir/2, and numerically we see that Ao — 3.68078. 

If 1 < A < Ao, then for t > we have ir < £1 < 2ir and < £2 — £1 < tt/2, and this implies that f(t) > for any 
t > 0. Similar inequalities show that f(t) > for t < 0, and we also have /(0) = 4. □ 

Remark. Of course, it is not necessary to have a function f{t) positive for all t, in order to obtain a positive integral. 
Numerically, one can see that the result of Proposition [21] is valid for other values of A, larger than Aq. 



Finally we point out, for integer values of A, the Mel'nikov potential (91 ) could be computed explicitly by writin 



it as the integral of a rational function. Applying standard trigonometric and hyperbolic formulas, and replacing 
t —> § + 1, the integral becomes 



L(s) 



sinh (I — t) 



sinh A (7 



t 



cosh (I — t) ■ cosh A (f + 1) 



dt. 



This can be transformed into the integral of a rational function, if A is integer, through the change x = tanht. This 
change is analogous to the one carried out : 
rational coefficients (in fact, the change x 



change is analogous to the one carried out in Section 3.2 transforming the Riccati equation (68) into an equation with 

tanht is the composition of q\ = 4arctane* and x = — cos(q\/2)). 



For instance, in the simplest case A = 1 we obtain 



L(s) 



! sinh 



2 s 



d.r 



'-1 (cosh 2 



: sinh 2 §)' 



= -4tanh| • ^ — 

2 I cosh 2 S 



2 tanh ■ 



A numerical inspection of this function (see Figure [2]) shows that it has 2 additional critical points, associated for 
e ^ small enough to other perturbed loops. 
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